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Abstract: We show that the algorithm originally proposed by Collins and Knowles 
for spin correlations in the QCD parton shower can be used in order to include spin 
correlations between the production and decay of heavy particles in Monte Carlo 
event generators. This allows correlations to be included while maintaining the 
step-by-step approach of the Monte Carlo event generation process. We present ex- 
amples of this approach for both the Standard and Minimal Supersymmetric Stan- 
dard Models. A merger of this algorithm and that used in the parton shower is 
discussed in order to include all correlations in the perturbative phase of event gen- 
eration. Finally we present all the results needed to implement this algorithm for 
the Standard and Minimal Supersymmetric Standard Models. 
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1. Introduction 

In modem particle physics experiments it is important to have a Monte Carlo simula- 
tion which accurately predicts the results of the experiment for both those processes 
which have already been observed and any new physics which may be discovered. 
As higher energies are probed this often means the production of heavy particles 
which decay before hadronizing giving decay products which are detected, e.g. the 
top quark. 

As these particles are often fermions the distributions of the decay products 
are affected by correlations between the production and decay of the fermion. In 
most Monte Carlo event generators these correlations are neglected. In the Standard 
Model (SM) this only applies to the top quark and tau lepton. However, in most 
models of new physics, for example supersymmetry, there are heavy fermions which 
decay and correlations between the production and decay of these fermions may be 
important [1-4]. 

There have been a number of calculations of the spin correlation effects for 
specific processes, but these calculations all require the hard collision process and the 
decay of the heavy fermions to be performed in the same step which is problematic 
in Monte Carlo event generators. 

There are however methods [5-8] which have previously been used in order to 
correlate all the partons produced in the QCD parton shower. In this paper we shall 
show how a very similar algorithm can be used in order to correlate the spins of all 
the heavy particle decays in the event while still maintaining both the step-by-step 
approach of the Monte Carlo event generator and an algorithm whose complexity 
grows only linearly with the number of particles. 

In the next section we will review the ideas of Monte Carlo simulations and 
discuss the current implementation of heavy particle decays in event generators. 
This is followed by a discussion of the algorithm we will use in Section ^ We 
provide several examples of the results of this algorithm for both Standard Model and 
supersymmetric (SUSY) processes in Section |[ We then discuss how to incorporate 
spin correlations in the decay of heavy particles as well as the spin correlations inside 
and between jets which are already included in some Monte Carlo event generators 
[9, 10] using the algorithm of [5-8]. Finally, we present the necessary results in 
order to implement this algorithm for both Standard and Minimal Supersymmetric 
Standard Model processes. 

2. Monte Carlo Event Generation 

There are a number of general purpose Monte Carlo event generators available [9- 
13]. The structure of the event generation procedure is basically the same in all 
these programs. The differences are in the algorithms used in the different steps of 
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generating the event. In general the Monte Carlo event generation process can be 
divided into five phases: 

1. The hard process where the particles in the hard collision and their momenta 
are generated, usually according to the leading-order matrix element. This 
can be of either the incoming fundamental particles in lepton collisions or of a 
parton extracted from a hadron in hadron-initiated processes. In the example 
event shown in Fig. ^ the hard process is e"^e~ — > tt. 

2. The parton-shower phase where the coloured particles in the event are per- 
turbatively evolved from the hard scale of the collision to the infrared cut-off. 
This is done for both the particles produced in the collision, the final-state 
shower, and the initial partons involved in the collision for processes with in- 
coming hadrons, the initial-state shower. In the example shown in Fig.|l| the 
top quarks radiate gluons and the gluons branch. 

3. Those particles which decay before hadronizing, e.g. the top quark and SUSY 
particles, are then decayed. Any coloured particles produced in these secondary 
decays are evolved by the parton-shower algorithm. These decays are usually 
performed according to a calculated branching ratio and often use a matrix 
element to give the momenta of the decay products. The example in Fig.|I| 
shows the semi-leptonic decay of both top quarks. 

4. A hadronization phase in which the partons left after the perturbative evo- 
lution are formed into the observed hadrons. For processes with hadrons in 
the initial state after the removal of the partons in the hard process, we are 
left with a hadron remnant. This remnant is also formed into hadrons by the 
hadronization model. In the example shown in Fig.|l] the cluster model [14], 
which is used in HERWIG [9, 10], is shown. 

5. Those unstable hadrons which are produced in the hadronization phase must 
also be decayed. These decays are usually performed using the experimentally 
measured branching ratios and a phase-space distribution for the momenta 
of the decay products. Any coloured particles produced in these decays are 
evolved according to the parton-shower algorithm and hadronized. This pro- 
cedure is repeated until all the unstable particles have been decayed. 

In the Standard Model the only quark which decays before hadronizing is the 
top quark. It is therefore possible to include tt production, including spin effects, as 
a 2 — > 6 process, i.e. including the decay of the top quarks, in a Monte Carlo event 
generator. However this leads to problems when including QCD radiation from the 
top quark before it decays. 

In all Monte Carlo event generators the particles are produced and then QCD 
radiation is generated via the parton-shower algorithm. Any unstable particles are 
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Figure 1: Example of a Monte Carlo event. This example shows the production of tt 
in e'''e~ collisions followed by the semi-leptonic decay of the top quarks. The cluster 
hadronization model [14] is shown where the gluons left after the parton-shower phase are 
non-perturbatively split into quark-antiquark pairs. The quarks and antiquarks are then 
paired into colour-singlet clusters using the colour flow information in the event. These 
clusters decay to give the observed hadrons. 

then decayed. This is the main problem for the inclusion of spin correlations, i.e. 
we want to be able generate the production and decay of the particles as separate 
processes rather than in one step so that QCD radiation can be generated for the 
particles before they decay. 

The problem is somewhat different for the inclusion of tau decays due to the large 
number of tau decay modes. This makes it impossible to include all the possible decay 
channels as 2 ^ n body processes. 

The situation in SUSY models is similar to that for tau decays. In order to 
facilitate the experimental search for supersymmetry a number of Monte Carlo event 
generators have either been extended to include supersymmetric processes [9-13,15, 
16], or written specifically for the study of supersymmetry [17,18]. 

These event generators have become increasingly sophisticated in their treatment 
of supersymmetric processes. In general, the decays of the supersymmetric particles 
produced in the hard collision process are assumed to take place independently. 
The event generators differ in how these decays are performed. While HERWIG 
[9,10,16] continues to use a phase-space distribution for the three-body decays of the 
supersymmetric particles ISAJET [13], PYTHIA [11,12,15] and SUSYGEN [17,18] 
use the full three-body matrix element. 

A number of studies have also been performed in which spin correlation effects are 
included, and these have been shown to be important for a future linear collider [1-4]. 
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Some spin correlation effects have been included in SUSYGEN [17, 18]. However in 
order to accomplished this both the production process and decays must be generated 
at the same time which limits both the number of processes and final-state particles 
that can be simulated in this way. 

In general SUSY models, particularly in hadron collisions, very complicated de- 
cay chains can occur and it is impossible to generate them all as 2 — ri particle 
processes. We therefore need an algorithm which both preserves the step-by-step 
approach of the traditional Monte Carlo event generators and a complexity which 
does not grow exponentially with the number of final-state particles. 

In this paper we will show that a completely general algorithm can be imple- 
mented which takes all these effects into account while still maintaining the step-by- 
step approach of the Monte Carlo event generation process. The results from the 
implementation of this algorithm in the HERWIG event generator are then compared 
to matrix element calculations from a number of observables. 

3. Spin Correlation Algorithm 

The spin correlation algorithm we will use is essentially identical to that presented 
in [5-8] for the case of spin correlations in the QCD parton-shower phase of the 
Monte Carlo event generation process. However, the algorithm presented in [5-8] 
only considers the case of 1 — 2 branchings, as these are all that occur in the QCD 
parton shower. 

The algorithm can be formulated entirely in terms of the matrix elements for the 
hard process and the various decays, spin density and decay matrices. The algorithm 
is defined as follows: 

1. The momenta of the particles in the hard process are generated according to 
the matrix element^ 

Pki<Pk24-^«i«:2;Ai...a„A1*;4;A'^...a;, H ^'>^^K' ^^'^^ 

i=l,n 

where A^ki,k2;Ai...a„ is the matrix element for the 2 — n body process, is 
the helicity of the ith incoming particle, Aj is the helicity of the ith outgoing 
particle, p^^; is the spin density matrix for the ith incoming particle and D\ y 
is the decay matrix for the ith outgoing particle. 

In the initial stage of the algorithm the momenta of the particles involved 
in the hard collision are generated according to Eqn.^Tl] with D\_y = Sx-y, 

^Throughout this paper we will use the Einstein summation convention where repeated indices 
are summed over. 
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and p* , = ^^k-k'. for unpolarized incoming particles. For polarized incoming 
particles 

_/|(l + P3) 

^^'^ V i(i-^3) 

where 7^3 is the component of the polarization vector parallel to the beam 
direction. For incoming antiparticles the sign of V3 must be changed. 

2. One of the outgoing particles is chosen at random and a spin density matrix 

P^i^'j " ^/'^i«'/«24-^'^i''2;Ai...A,...A„>t;^;4;A;...A'...A; JI ^A^AJ' ^^■^'> 

calculated for the decay of this particle. The normalization 

Np = P;li«'/«24-^«i«2;Ai...A,...A„At:;4;A;...A,...A;, H (3-4) 

is chosen so that the trace of the spin density matrix is one. 

3. The decay mode of this particle is selected according to the branching ratios 
and the momenta of the particles produced in the n-body decay generated 
according to 

Paoa;,A^Ao;Ai...a„M1-a;...a; H (3-5) 

i=l,n 

where Aq is the helicity of the decaying particle and Aj is the liclicity of the ith 
decay product. As before the decay matrices are taken to be D^XiX' ~ 
the initial step of the algorithm. 

4. One of the particles produced in this decay is selected and a spin density matrix 
calculated. Again the normalization 

Ndp = PAoAJ,A^Ao;Ai...A,...A„A^I-a',..A,...A; n ^A,A^' (3-7) 

is chosen such that the spin density matrix has unit trace. This spin density 
matrix is used as the input to the third step of the algorithm. 

5. The third and fourth steps are repeated until the particle selected in step four 
is stable. When this occurs the decay matrix for this stable particle is set to 
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D\,y = 6x.x'^ and another particle from that decay selected to be decayed next 
using a spin density matrix for this particle calculated using Eqn.|3]^. When 
all the particles produced in a given decay have been developed a decay matrix 

^AoA^ = ^^Ao;A....A„A<I;,V^..... H ^^A^, (3-8) 

i=l,n 

is calculated for the decay. Again the normalization 

Nn = -MAo;A....A„-MI„;Ai...A; H ^A.A^' (^.O) 

i=l,n 

is chosen such that the trace of the decay matrix is one. 

6. A new particle is selected from the decay which produced the decaying particle. 
This new particle has a spin density matrix given by Eqn.^]^ using the decay 
matrix calculated from Eqn. |3.8| for the particles which have already been de- 
cayed rather than the identity. In this way the decay products of this particle 
will have the correct correlations with the decay products of the other parti- 
cles produced in the same decay. This step is repeated until all the particles 
produced in a decay are developed. Eqn. ^]8| is then used to compute the decay 
matrix for this decay and the previous decay in the chain is developed. 

7. Eventually this will give a decay matrix for the particle produced in the hard 
process. At this point a new particle from the hard process is selected to 
be decayed with a spin density matrix given by Eqn.^]^ with the identity 
replaced by the calculated decay matrix for those particles which have already 
been decayed. In this way the correlations between the decays of the particles 
produced in the hard process are generated. This is repeated until all the 
particles produced in the hard process have been decayed. 

At each point in the algorithm the normalization is chosen such that the trace 
of the spin density and decay matrices is one. This is necessary in order to maintain 
the probabilistic interpretation of the spin density matrices and proves to be usefully 
for the decay matrices. 

4. Examples 

There are many quantities for which it is interesting to calculate the spin correlations. 
In this section we will merely show a few examples of the method we propose for 
including the spin correlations and comparisons with the full analytic result for these 
quantities. Due to the complexities in calculating these full results this will limit us to 
relatively simple observables, although in principle much more complicated quantities 
can be calculated. The quantities have been chosen to both illustrate the technique 
and be phenomenogically important. 
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Figure 2: Feynman Diagrams for the decay qj, qX2 ~^ Q^^^/j' 
4.1 ^ qX° ^ q^±^l 

In hadron-hadron collisions the coloured SUSY particles, i. e. the squarks and gluinos, 
are preferentially produced and therefore the main source of electroweak gaugino and 
slepton production can be in the cascade decays of the coloured sparticles. 

In many models the decay chain — > qxl ^ l^^^i? ~^ Fig-S is 

important because by studying edges in the £~^£~ , and qi~^£~ mass distributions 
the masses of the lightest two neutralinos, squark and slepton can be reconstructed 
[19-23]. 

However, if we consider the spins of the particles involved the quark produced 
in the squark decay will be left-handed in the massless limit. The lepton produced 
in the neutralino decay will be right-handed. This will lead to very different decay 
distributions for the charge conjugate decay modes of the neutralino. 

This difference is most noticeable in distribution of the mass of the quark and the 
lepton produced in the neutralino decay. Fig. |. The results in Fig. || were generated 
at SUGRA point 5 [20,24], i.e. Mq = 200 GeV, M1/2 = 300 GeV, Aq = 300 GeV, 
tan (3 = 2.1 and sgn where Mq is the universal SUSY breaking scalar mass at the 
Grand Unified Theory (GUT) scale, M1/2 is the universal SUSY breaking gaugino 
mass at the GUT scale, Aq is the universal tri-linear soft SUSY breaking parameter 
at the GUT scale, sgn/x is the sign of the parameter and tan/5 is the ratio of 
the two Higgs vacuum expectation values. At this SUGRA point the squark mass 
Mq^ = 833.8 GeV, the lightest neutralino mass M^o = 115.5 GeV and the next- 
to-lightest neutralino mass M^o = 213.2 GeV. The SUSY spectrum was generated 
using ISAJET7.51 [13]. 

The differences in the shapes of the two distributions in Fig. | can be understood 
by considering the helicities of the particles, the quark is left-handed while the pro- 
duced lepton/antilepton is right-handed. In the case of an antilepton this means that 
if the antilepton and the quark are back-to-back they have no net spin and therefore 
this configuration which gives the edge in the mass distribution shown in Fig.^ is 
favoured. However, if a lepton is produced this configuration has net spin one and 
therefore cannot be produced in the decay of a scalar so the distribution in Fig. ^b 
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Figure 3: Distribution of the mass of quark and lepton produced in the decay 
~^ qX2 q-^^^ij- The sohd hne gives the result of phase space, the dashed hne gives 
the full result and the dotted line the result of the spin correlation algorithm. 

vanishes at the end point. 

Unfortunately in any experiment there is no way to distinguish between a quark 
and an antiquark. Hence the average of the distributions in Fig. |^ will be observed. 
This average is in good agreement with the phase-space distribution. 

This is the simplest, non-trivial, application of the spin correlation algorithm. 
Here all that is necessary is that the spin density matrix for the production of the 
neutralino in the scalar decay is used to perform the decay of the neutralino to a 
scalar and a fermion. The matrix elements are therefore relatively simple. Indeed 
the matrix elements are simple enough to allow us to compare the results of the 
algorithm and the full result analytically. The matrix element for the 3-body process 
Ql qX2 ^^'^^R is given by 

-Mq,^q^o^q,+|- = u{p^)a'Px {A - ^x) h^'PMvd. (4.1) 

where is the four-momentum of the neutralino, Pq is the four-momentum of the 
quark, is the four-momentum of the produced lepton. Pa = |(1 + A75) , and 
are the couplings for the production and decay of the neutralino, respectively. In the 
notation used in Appendix El, = a*~^- and h'^ = a\.r „ where in the limit that 
we neglect left/right sfermion mixing = 5a,p=i for left squarks and L^^ = 5a,(3=2 
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for right sleptons. As the width of the neutrahno is small compared to its mass we 
have assumed that it is on-mass-shell. 

We can now consider how the spin correlation algorithm attempts to reproduce 
this result. The matrix element for the first step of the process is 



M 



(4.2a) 



where Xy. is the helicity of the neutralino, Aq is the helicity of the quark, Zq and 
are the reference vectors for the quark and neutralino, respectively. The use of 
these reference vectors in defining the spinors for massive fermions is described in 
Appendix 0. The sign of the helicity of the neutralino has be chosen in order for the 
outgoing neutralino to be considered as a particle rather than an antiparticle. 

In the first stage of the algorithm the momenta of the neutralino and quark are 
generated according to the spin averaged matrix element 



1+2 _2 



2a^a mc^niy^ 



(4.3) 



spins 



As the quark is stable we can average over its spin and produce the spin density 
matrix needed to perform the decay of the neutralino, 

1 1 



^qO"^ + rriqa^a^ 



(4.4) 



where the normalization is chosen to give YlxP^>' ~ ^ ^^"-^ therefore equal to 
the spin averaged matrix element, Eqn. [4.3| . 

The matrix element for the decay of the neutralino is given by 



X2- 



1 



--ux^ik) (/x -"^x) 



(4.5a) 

mib-^']ux,ik), (4.5b) 



^VPx ■ kPi ■ 

where p^ is the four-momentum of the lepton, is the reference vector used to define 
the direction of the lepton's spin and Xi is the helicity of the lepton. 

As with the quark because the produced lepton is stable we can average over its 
helicities giving 

1 



2p^ ■ 



-ux^ik) iA-^x) 



(4.6) 



iA-mx)ux'(lx)- 
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I I 
e+ < ^ ► X? e+ < ^ ► 

Figure 4: Feynman diagrams for the production of X2X1 iii e~'^e~ collisions. 

This can now be contracted with the spin density matrix and the sum over the 
hehcities of the neutrahno performed giving 

— Tr |(/^ - m^) ^/qfl"^ + mqO+a" j P+ + ^Ki^"^^ + mqa+a" j P_ 
(^^ - m^) (^^^6+^ - m^6+r ) P_ + (^/^r ^ - m^6+r ) P+ } . 

The second decay is generated according to this formula. The normahzation of the 
spin density matrix is cancelled by the matrix element which is used to generate 
the first decay. If we compare the above result with the matrix element for the full 
3-body process we can see that it agrees with the amplitude squared. Hence, the 
decay products will have the same distribution as the full three-body matrix element. 
This can be seen in Fig. ^ where the result of the full 3-body matrix element and the 
spin correlation algorithm are in very good agreement. 

4.2 e+e- ^ xUi 

In the MSSM the lightest supersymmetric particle (LSP), usually taken to be the 
lightest neutralino, is stable and weakly interacting. It therefore escapes from the 
detector without interacting giving missing transverse energy. This means that in 
a future linear collider the production of X2X1) Fig-@) ^lay well be the detectable 
supersymmetric final state which requires the smallest centre-of-mass energy. The 
production of the second-to-lightest-neutralino will be followed by its decay to the 
lightest neutralino. This decay will be X2 ~^ Xi^f via either real or virtual sfermions, 
Higgs or Z bosons, Fig.|^. 

There have been a number of studies of spin correlations for this process both 
with [3,4] and without [1] beam polarization. The study of both the decay correla- 
tions and polarization effects is important because it would allow the nature of the 
neutralinos to be determined. 

This is a more complicated example of the use of the spin correlation algorithm 
than the squark decay we studied in the previous section. Instead of being produced 
in a scalar decay the neutralino is produced in a 2 ^ 2 process where there are a 
number of Feynman diagrams and the possibility of polarizing the incoming particles. 
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Figure 5: Feynman diagrams for the decay X2 ~^ ffXi- The exchange of the MSSM Higgs 
bosons is only important for the third generation of fermions. 



However, as before the spin correlation algorithm only has to compute the spin 
density matrix for the production and use this to generate the decay of the neutralino 
which may be via either a two or three body process, Fig. ^. 

The easiest case to consider first is the decay of the neutralino to a lepton and 
a right-handed slepton, as in Section [4.1| , again using SUGRA Point 5. As can be 
seen in Fig. ^ there is a correlation between the direction of the produced lepton and 
the beam direction and this correlation is affected by the polarization of the beam. 
It should be noted that fully polarized beams are not achievable in practice and are 
only shown in order to illustrate the results of the algorithm. As with the results in 
Section the result of the spin correlation algorithm is in good agreement with the 
result from a 4-body matrix element including the decay of both the neutralino and 
slepton. 

It is also possible that neither the gauge boson or sleptons in the diagrams in 
Fig.^ can be real. As an example we considered a SUGRA point with non-universal 
gaugino masses at the GUT scale in order to decrease the mass difference between 
the lightest two neutralinos. We used the point Mq = 210 GeV, Aq = GeV, 
tan p = 10, Ml = 450 GeV, M2 = 350 GeV, M3 = 350 GeV, where Mi is the soft 
SUSY breaking mass for the bino at the GUT scale, M2 is the soft SUSY breaking 
mass for the wino at the GUT scale and M3 is the soft SUSY breaking mass for the 
gluino at the GUT scale. At this point the lightest neutralino is dominantly bino-like 
and the next-to-lightest neutralino is dominantly wino-like, the lightest neutralino 
mass is M^o = 182.4 GeV, the next-to-lightest neutralino mass is M^o = 264.2 GeV, 
the right-slepton mass is M^^ = 271.3 GeV and the left-slepton mass is M^^ = 
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(a) Unpolarized 



(b) 6,6] 



(b) 6j;6^ 






Figure 6: Angle between the lepton produced in e'''e~ X2X1 ~^ ^r^~Xi the incom- 
ing electron beam in the laboratory frame for a centre-of-mass energy of 500 GeV with (a) 
no polarization, (b) negatively polarized electrons and positively polarized positrons and 
(c) positively polarized electrons and negatively polarized positrons. The solid line shows 
the default result from HERWIG which treats the production and decay as independent, 
the dashed line gives the full result from the 4-body matrix element and the dotted line 
the result of the spin correlation algorithm. 



322.5 GeV. This point was chosen so that both the right-handed slepton and the Z 
boson could be almost on mass-shell in the three body decay X2 ~^ Xi^~^^~ ■ 

The angle between the lepton produced in the X2 decay and the beam is shown in 
Fig.|^ for a three choices of polarization for the incoming beams. The angle between 
the produced leptons is shown in Fig. ^ for the same choices of beam polarization. 
The default treatment of this process in the HERWIG event generator is that the 
production and decay of the next-to-lightest neutralino take place independently and 
the neutrahno is decayed using a phase-space distribution for the decay products. 
We have therefore shown the result of using a matrix element to perform this decay, 
while still treating the production and decay as independent, as well as the result of 
the spin correlation algorithm and a full calculation of the process in Figs. ^ and p. 

As can be seen in Fig.|^ the inclusion of the matrix element for the decay on 
its own has no effect on the angle between the lepton and the beam, whereas the 
inclusion of this matrix element significantly improves the agreement between the 
result of HERWIG and the full result for the distribution of the angles between the 
produced lepton and antilepton. Fig. |[ Indeed the inclusion of the matrix element 
seems to be all that is required in order for the Monte Carlo simulation to reproduce 
the full result for this distribution. The full spin correlation algorithm is necessary 
to reproduce the correlation between the beam direction and the direction of the 
produced lepton. Again, there is good agreement between the full 4-body matrix 
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(a) Unpolarized 
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Figure 7: Angle between the lepton produced in e"'"e~ — > X2X1 ~^ '^~^^~XiXi and the in- 
coming electron beam in the laboratory frame for a centre-of-mass energy of 500 GeV with 
(a) no polarization, (b) negatively polarized electrons and positively polarized positrons 
and (c) positively polarized electrons and negatively polarized positrons. The solid line 
shows the default result from HERWIG which treats the production and decay as inde- 
pendent and uses a phase-space distribution for the decay products of the neutralino, the 
dot-dashed line also includes a matrix element for the neutralino decay, the dashed line 
gives the full result from the 4-body matrix element and the dotted line the result of the 
spin correlation algorithm. 




Figure 8: Angle between the lepton and antilepton produced in 
e''~e~ X2X1 ~^ ^~^^~XiXi in the laboratory frame for a centre-of-mass energy of 
500 GeV with (a) no polarization, (b) negatively polarized electrons and positively polar- 
ized positrons and (c) positively polarized electrons and negatively polarized positrons. 
The lines are described in the caption of Fig.|^. 
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Figure 9: Angle between the lepton produced in e''~e~ ^ tt — > hhi~^ i^e and the incom- 
ing electron beam in the laboratory frame for a centre-of-mass energy of 500 GeV with (a) 
no polarization, (b) negatively polarized electrons and positively polarized positrons and 
(c) positively polarized electrons and negatively polarized positrons. The solid line shows 
the default result from HERWIG which treats the production and decay as independent 
but includes a matrix element for the weak decay of the top, the dashed line gives the full 
result from the 6-body matrix element, the dotted line the result of the spin correlation 
algorithm and the dot-dashed line the result of the spin correlation algorithm when the 
decay matrix for the first quark is neglected. 



element result and the spin correlation algorithm for both distributions. 
4.3 e+e- tt 

A future linear collider will have both the energy to produce top quark pairs and the 
possibility of polarized incoming beams. An accurate measurement of the top quark 
mass by either scanning the threshold or reconstructing the decay products will be a 
goal of such an experiment. There have been a number of studies of spin correlations 
in top quark production at such a machine [25]. 

Here we are mainly concerned with this process as an example of the application 
of the spin correlation algorithm. In this process in addition to the complications of 
the 2 — i> 2 process, which can have polarization of the incoming particles, there is the 
problem of correlating the decay of the top quarks with each other. In addition to 
using the spin density matrix for the top production to perform the decays the decay 
matrix from the first decay must be used in calculating the spin density matrix used 
for the decay of the second quark. In the default treatment of this process HERWIG 
performs the production and decay of the quarks independently but uses the full 
3-body matrix element for the quark decays. 

Fig.|^ shows the angle of the lepton produced in the top antiquark decay with 
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(a) Unpolarized 





(b) e^el 




Figure 10: Angle between the lepton and top quark produced in e'''e~ — > tt — > hhi~^ i'£i~ 
in the laboratory frame for a centre-of-mass energy of 500 GeV with (a) no polarization, 
(b) negatively polarized electrons and positively polarized positrons and (c) positively 
polarized electrons and negatively polarized positrons. The lines are described in the 
caption of Fig. ^. 




Figure 11: Angle between the leptons produced in e''~e~ ^ tt — > hhl~^ i'£l~ i>£ in the labora- 
tory frame for a centre-of-mass energy of 500 GeV with (a) no polarization, (b) negatively 
polarized electrons and positively polarized positrons and (c) positively polarized electrons 
and negatively polarized positrons. The lines are described in the caption of Fig. ^. 



respect to the incoming electron beam. Fig. |T0| shows the angle between the lepton 
and the top quark and Fig. ^ shows the angle between the lepton and the antilepton. 
The result of the spin correlation algorithm agrees well with the full 6-body matrix 
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element for all these distributions. Again the major discrepancies between the default 
treatment in HERWIG and the full result are for the correlation of the direction of 
the outgoing leptons with the beam and the correlation between the direction of the 
produced lepton and the second quark in the event. 

In addition to the full result of the spin correlation algorithm we have also 
included the results of the algorithm where the identity is used rather than the decay 
matrix calculated for the decay of the first quark when calculating the spin density 
matrix for the decay of the second quark. Neglecting this step of the algorithm 
still gives good agreement with the matrix-element calculation for the correlation 
between the direction of the produced lepton and the beam direction, Fig. |], and 
the correlation between the direction of the top quark and the lepton, Fig. This 
is because these correlations mainly depend on the details of the hard production 
process and not how the first quark decayed. 

However, in order to reproduce the result of the matrix element for the angle 
between the produced lepton and antilepton. Fig. 0, it is essential to include the 
decay matrix from the first quark as this is a correlation between the two decays. 
This effect is most noticeable when there is no beam polarization and both possible 
initial-state spin configurations contribute to the result. 

4.4 tt production in Hadron-Hadron Collisions 

Since the discovery of the top quark at the Tevatron there have been a number 
of studies of spin correlations in top quark pair production for both Run II of the 
Tevatron and the LHC There are two processes which contribute to tt production 
in hadron collisions 

qq — > tt, (4.8a) 
gg tt. (4.8b) 

The qq annihilation process presents no additional complications and is very similar 
to the e"'"e~ annihilation we discussed in the previous section. 

However the colour flows in the three diagrams for tt production via gluon scat- 
tering, Fig.|T^, are very different. It is easiest to extract the colour matrices from 
the helicity amplitudes and perform the colour sum/averages separately. The colour 



matrices for the three diagrams, shown in Fig. 0, are 



CtcL = tlA,^^ t-channel, (4.9a) 



Cute, = t'c^c'K'c, n-channel, (4.9b) 

Cstc, = ^r'X.c, s-channel, (4.9c) 

where /"^^ is the SU{3) structure constant, is the SU{3) colour matrix in the 
fundamental representation, Cs^t,u'tc2 colour matrices for the s-, t- and u- 

channel colour flows respectively, a is the colour of the first incoming gluon, b is the 
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Figure 12: Feynman diagrams for gg — > tt. 

colour of the second incoming gluon, ci is the colour of the outgoing quark, d is the 
colour of the intermediate quark and is the colour of the outgoing antiquark. 

In general the colour flow for diagrams involving the triple gluon vertex, i.e. 
the s-channel diagram, is not unique and can be rewritten in terms of the t- and 
It-channel colour flows using the definition of the structure constants 

%f°^^t'c^C2 ~ ^'ciC'^c'C2 ~ ^C1C'^C'C2' (4-10) 

giving 

/~i ab ah ab ( A ^ ^\ 

SO that we only have two different colour flows to deal with. 

In principle the presence of different colour flows could be a significant problem 
to the algorithm we are using. We should keep track of the colour flow in the 
production processes and decays and calculate spin density matrices for each of the 
flows and sum over them. If the decays have one than one colour flow, in addition 
to the different colour flows in the hard process, this could lead to an increase in 
the complexity of the algorithm with the number of final state particles, i.e. the 
number of colour flows would grow exponentially. However, in both the Standard 
and Minimal Supersymmetric Standard Models, provided we only consider three- 
body decays as we are doing, there can only be different colour flows in the hard 
production process and it is therefore sufficient to replace Eqn. |3.1| with 

^ <^abP«i«'/^24-M°^^2.^^...^„A<t|^.V^ JJ L'^^y^, (4.12) 

a,b=l,C]\f i=l,n 

where A^*^ is the matrix element for colour flow a. Cat is the number of colour flows 
and Cab is the colour factor for the matrix elements. Eqns. p.3| and p.4| must also be 
modified in the same way. 
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For example in tt collisions 



1 



Ctt — (Ar2_i)2^cic'^c'c2^c"ci^C2c" " iiV"' (4.13a) 



c 



Cuu = (]v^3ip'^cic'^c'c2^c"ci^C2c" ~ iiV"' (4.13b) 

Ct« = (Nl-i)iKic'tc'c2^c"citc2c" = ~4jV'c(A^2 _ ly (4.13c) 

the colour factor Cut = Cut-"^ 

Another option which would be necessary if the decays have more than one 
possible colour flow, in for example R-parity violating SUSY models [26,27], is to 
select one of the colour flows in the same way as is already done for the production of 
QCD radiation [28]. In this procedure we would first select a colour flow, a, using [28] 

CaaP«,«iP',4A^:,.2;Ai...A„A^<.-Ai...A„^a, (4.14) 

where 

_ I^b,c=l,CAr ^''cP^jk;PLk2"^'^i''2;Ai...A„-^k;4;Ai...A„ (A'\r:\ 
l^b=l,CN ^^^-^ 'kiK2;Ai...A,/^'44;Ai...A„ 

This procedure ensures that the sum over the redefined colour flows reproduces the 
full result, Eqn. [4.12| . This is used to select the colour flow a for the process and 



instead of Eqn.^TT] to give the distributions of the particles produced in the hard 
process. The matrix element for the selected colour flow is then used with Eqn. ^]3| 
to calculate the spin density matrices for the decays of the particles produced in the 
hard process. 

In models where there is more than one colour flow in a decay we would select the 
colour flow for the decay and generate the momenta of the decay products according 
to 

PA„AJ,aaA^A,;A,...A„-^AU....A„^., (4-16) 

where 

p _ X]fe,c=l,Cjv '-'''c-^Ao;Ai...A„-^r^;Ai...A„ , - 

2^6=l,CAr '-^fc''./\/lAo;Ai...A„-/^'Aj,;Ai...A„ 

The matrix element for the selected colour flow is then used to calculate the spin 
density matrices for the particles produced in the decay and the decay matrix for 
this decay once all the particles produced in the decay have been developed. 



^The sign of Ct„ is different to that given in Table|T^ because it proved easier to absorb this sign 
as a change of sign of the u-channel matrix element. 
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01% 


01% 


Full 


Spin 1 


Spin 2 


Full 


Spin 1 


Spin 2 




2.46 ± 0.04 


2.45 ± 0.02 


2.41 ± 0.02 


2.42 ± 0.04 


2.46 ± 0.02 


2.44 ± 0.02 


02 


2.94 ±0.02 


2.92 ±0.01 


2.90 ±0.01 


2.91 ±0.02 


2.91 ±0.01 


2.89 ±0.01 


03 


1.68 ±0.02 


1.67 ±0.01 


1.64 ±0.01 


1.67 ±0.02 


1.68 ±0.01 


1.65 ±0.01 


04 


0.03 ±0.02 


0.01 ±0.01 


0.01 ±0.01 


-0.01 ±0.02 


0.00 ±0.01 


-0.01 ±0.01 


05 


1.95 ±0.02 


1.94 ±0.01 


1.92 ±0.01 


1.92 ±0.02 


1.93 ±0.01 


1.92 ± 0.01 


0, 


-0.48 ±0.03 


-0.47 ±0.02 


-0.46 ±0.02 


-0.49 ±0.03 


-0.45 ±0.02 


-0.48 ±0.02 



Table 1: Average values of the observables 0\-% and (5i_6 at Run II of the Tevatron. 
The top (antitop) quark was required to have rapidity |y| < 2.0 and transverse momentum 
■pT > 15GeV for the observables Oi-q and Oi-q, respectively. Spin 1 gives the results 
of the spin correlation algorithm using Eqn. 4.12 and spin 2 gives the results of the spin 



correlation algorithm using Eqn. 4.14 



This approach neglects the interference between the different colour flows in 
generating the spin correlations, just as they are neglected when generating the QCD 
radiation. As this procedure is not necessary for the processes we are considering we 
will use the first procedure in processes where the hard process has more than one 
colour flow. 

There have been many studies of spin correlations in hadron-hadron collisions. 
In order to study the results of the spin correlation algorithm we have chosen to 
follow the approach of [29,30] in which the following observables were defined for top 



quark pair production followed by t ^ bW"^ bqq and t —>■ bW 



hi-j?: 



(Qb ■ Pp) (q^- ■ Vp] 

qfe ■ kt) ■ kt ; , 

1 - 

2 . 

(pp X qD ■ (pp X q^-) = Ci - C2 



qfe ■ Pp) q£- ■ kt + q^ ■ ki ) (q^- ■ Pp 



(4.18a) 
(4.18b) 
(4.18c) 

(4.18d) 
(4.18e) 
(4.18f) 



All the hatted vectors are unit vectors and the quantities without an asterisk are 
measured in the laboratory frame, Pp is the beam direction, q^- is the direction of 
the charged lepton, is the direction of the top quark and yt is the rapidity of the 
top quark. The asterisked quantities are measured in the top quark rest frame and q^ 
is the direction of the bottom quark. In a similar way the quantities (5i_6 are defined 
for the charge conjugate decay modes. The observable O^^ is zero when integrated 
over a symmetric rapidity interval which is why the observable is defined. 

The results for these observables at Run II of the Tevatron with a centre-of- 
mass energy of 2TeV are shown in Table I] where for events with t bW"^ bqq 
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0/% 


0/% 


Full 


Spin 1 


Spin 2 


Full 


Spin 1 


Spin 2 




-2.01 ± 0.02 


-2.04 ± 0.02 


-2.04 ± 0.02 


-2.03 ± 0.02 


-2.04 ± 0.02 


-2.02 ± 0.02 


02 


-0.16 ±0.02 


-0.16 ±0.02 


-0.18 ±0.02 


-0.18 ±0.02 


-0.19 ±0.02 


-0.16 ±0.02 


03 


-0.55 ±0.02 


-0.56 ±0.02 


-0.57 ±0.01 


-0.57 ±0.02 


-0.59 ±0.02 


-0.57 ±0.01 


04 


0.00 ±0.02 


-0.01 ±0.01 


-0.02 ±0.01 


-0.02 ±0.02 


0.00 ±0.01 


0.01 ±0.01 


05 


-0.23 ±0.02 


-0.25 ±0.01 


-0.25 ±0.01 


-0.26 ±0.02 


-0.28 ±0.01 


-0.25 ±0.01 




-1.85 ±0.02 


-1.87 ±0.02 


-1.86 ±0.01 


-1.85 ±0.02 


-1.85 ±0.02 


-1.85 ±0.01 



Table 2: Average values of the observables Oi-q and Oi-q at the LHC. The top (antitop) 
quark was required to have rapidity |y| < 3.0 and transverse momentum pT > 20 GeV for 
the observables Oi_6 and Oi-q, respectively. Spin 1 gives the results of the spin correlation 
algorithm using Eqn. 4.12 and spin 2 gives the results of the spin correlation algorithm using 
Eqn.|4.14. 



and t —>■ bW~ —>■ hi~iy we have required the top quark to have rapidity \yt\ < 2.0 
and transverse momentum pt > 15 GeV. For the charge conjugate decay modes the 
antitop quark was required to have rapidity \yt\ < 2.0 and transverse momentum 
Pt > 15 GeV. 

The results for the same observables at the LHC with a centre-of-mass energy of 
14 TeV are shown in Table 0. For events with t — > bW^ bqq and t bW~ hi~u 
we have required the top quark to have rapidity \yt\ < 3.0 and transverse momentum 
Pt > 20 GeV. For the charge conjugate decay modes the antitop quark was required 
to have rapidity \yt\ < 3.0 and transverse momentum pt > 20 GeV. 

All the results were generated using a top quark mass of 175 GeV and the default 
parton distribution functions in HERWIG6.3 [10] which are the average of the central 
and higher gluon leading-order fits of [31]. 

As can be seen in both Tables|l| and ^ the results of the spin correlation algorithm, 
with either method of handling the different colour flows, are in good agreement with 
the result of the 6-body matrix element. The default result from HERWIG for all 
these observables was consistent with zero. This is because these observables are 
designed to be sensitive to the spin correlations and zero in the absence of such 
correlations. 

At the Tevatron the most significant of these observables is O2 whereas the 
observable Oq is the most important at the LHC. At the LHC this is because while 
the average value of Oi is greater the fluctuations are also greater which gives a 
greater statistical sensitivity for Oq [29,30]. 

4.5 e+e- 

So far we have only considered processes in which the fermions have small widths 
relative to their masses, the next-to-lightest neutralino in most SUSY models has a 
width of less than one hundred MeV and the top quark width is about one GeV. We 
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Figure 13: Angle between the lepton produced in e'''e~ — > Xi'xF ~^ i'eXi^~ '^^Xi 
incoming electron beam in the laboratory frame for a centre-of-mass energy of 500 GeV with 
(a) no polarization, (b) negatively polarized electrons and positively polarized positrons 
and (c) positively polarized electrons and negatively polarized positrons. The solid line 
shows the default result from HERWIG, the dashed line gives the full result from the 
6-body matrix element and the dotted line the result of the spin correlation algorithm. 



(a) Unpolarized 



(a) Unpolarized 



(c) e^et 



0.4 — 






Figure 14: Angle between the leptons produced in e''~e^ XiXi ~^ i^eXi^'^eXi in 
the laboratory frame for a centre-of-mass energy of 500 GeV with (a) no polarization, 
(b) negatively polarized electrons and positively polarized positrons and (c) positively 
polarized electrons and negatively polarized positrons. The solid line shows the default 
result from HERWIG, the dashed line gives the full result from the 6-body matrix element 
and the dotted line the result of the spin correlation algorithm. 
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can therefore neglect the width in top quark and neutrahno production processes. In 
the hmit that the width of the decaying fermion can be neglected the helicity ampli- 
tudes for the production processes, including the fermion decays, can be factorized 
into a production and decay helicity amplitude. The spin correlation algorithm is in 
good agreement with the full matrix element in these cases. 

However, if the decaying fermion is not on mass-shell the process cannot be 
factorized and the result of the spin correlation algorithm may not give such good 
agreement with the 2 — >■ n body matrix elements. This should not be a serious 
problem because even in SUSY models the electroweak gauginos tend to have widths 
of at most a few GeV, for the heavier gauginos. The gluino width tends to be larger, 
but the gluino is heavier, and therefore the effect is the same. The particles which 
usually have the largest widths in most SUSY models are the squarks and these 
effects can be taken into account because they are scalars. 

In order to study the effects of the widths of the decaying fermions it is easiest 
to study chargino production in e"'"e~ collisions because this allows us to consider 
the production of xfx^ in which the lighter chargino has a very small width and the 
heavier chargino has a much greater width. Spin and polarization effects in chargino 
pair production in e"'"e~ collisions have been previously considered in [2,4]. 

We used the following SUSY parameters in this study: Mi = 78 GeV, 
Ma = 158 GeV, fi = 330 GeV, tan/5 = 3, M^^ = 171 GeV and M^^ = 180 GeV, 
where the soft SUSY breaking masses for the gauginos and the parameter are 
given at the electroweak scale, M/,^ is the soft SUSY breaking mass for the lepton 
doublet at the electroweak scale and Mg^ is the soft SUSY breaking mass for the 
right-handed electron singlet at the electroweak scale. 

At this point the lightest neutralino mass is M^o = 71.6 GeV, the lightest 
chargino mass is M-+ = 127.8 GeV, the heaviest chargino mass is M-+ = 357.6 GeV, 
the sneutrino mass is M^;^ = 160.9 GeV, the left selectron mass is M^^ = 176.1 GeV 
and the right selectron mass is M^^ = 184.2 GeV. The width of the lightest chargino 
is r -+ = 73 MeV and the heaviest chargino is T-+ = 2.9 GeV. 

This point was chosen to have a negligible lightest chargino width and a variety 
of leptonic decay modes, i.e. ~^ ^^^l ^~^^~Xi^ xt ~^ ^t^ ~^ ^~^^~Xi and 
xt Xi^^ ^^^~Xii of the heaviest chargino. 

It is interesting to first consider XiXi production, where there are no width 
effects. Fig. shows the angle between the lepton produced in e+e~ XiXi ~^ 
X^i^i>Xii~i^ and the beam direction and Fig. |l^ shows the angle between the lepton 
and the antilepton. As before there is good agreement between the full result and 
the spin correlation algorithm. 

In a Monte Carlo simulation however we need to include the widths of the un- 
stable particles as these widths can potentially be measured. When we include the 
effect of the width in the spin correlation algorithm we have two choices. We can 
use the off-mass shell momentum in Eqn. [A.2| for the spinor of the massive particle 
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Figure 15: Angle between the lepton produced in e'''e~ — > X1X2 ~^ i'eXi^~ '^^Xi 
incoming electron beam in the laboratory frame for a centre-of-mass energy of 500 GeV with 
(a) no polarization, (b) negatively polarized electrons and positively polarized positrons 
and (c) positively polarized electrons and negatively polarized positrons. The solid line 
shows the default result from HERWIG, the dashed line gives the full result from the 
6-body matrix element and the dotted line the result of the spin correlation algorithm. 



together with either the physical mass of the particle or the off-shell mass of the 
particle. 

In the first case this gives 

y 

^ma(p)ma(p) =#+m+ ^^(m^ (4.19) 

A 

where p is the off-shell four-momentum, m is the on-shell mass and / is the reference 
vector. This form reduces to the correct spin sum if the particle is on-mass shell. 
However for off-mass shell particles there is an additional component which depends 
on the choice of the reference vector and is related to how off-mass shell the particle 
is. 

If we choose to use the off-shell mass rather than the physical mass in Eqn. A. 2 
we obtain 

J2 ux{p)ux{p) = /+"^ofr, (4.20) 

A 

where rriog is the off-shell mass. Here the result again fails to reproduce the correct 
result by an amount which depends on the how far off-shell the particle is, however 
it does not depend on the arbitrary reference vector. We have studied the effects of 
both these choices. 
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Figure 16: Angle between the lepton produced in e'''e~ — > X1X2 ~^ i'eXi^~ '^^Xi 
incoming electron beam in the laboratory frame for a centre-of-mass energy of 500 GeV with 
(a) no polarization, (b) negatively polarized electrons and positively polarized positrons 
and (c) positively polarized electrons and negatively polarized positrons. The physical X2 
width was used. The solid line shows the full result, the dotted line shows the result of the 
spin correlation algorithm with the physical mass and the dashed line shows the result of 
the spin correlation algorithm with the off-shell mass of the chargino. 



(a) Unpolarized 





(c) e^et 




Figure 17 : Angle between the lepton produced in e''~e~ — > X1X2 ~^ ^eXi^~ ^£Xi ^^'^ the 
incoming electron beam in the laboratory frame for a centre-of-mass energy of 500 GeV with 
(a) no polarization, (b) negatively polarized electrons and positively polarized positrons 
and (c) positively polarized electrons and negatively polarized positrons. Ten times the 
physical X2 width was used. The lines are described in the caption on Fig. 
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We considered X1X2 production followed by the decay of the lightest chargino 
xf — ^ ^^z^Xi via a three-body decay and heaviest chargino xf ~^ W^Xi ~^ ^^^X^- 
The angle of the lepton produced in the X2 decay with respect to the beam direction 
is shown in Fig. |15| where we have neglected the widths of the decaying charginos. 

We then considered the effect of including the width of the decay chargino. 
Fig.|l^ shows the result for the physical width of the heaviest chargino and Fig.|l7 
shows the effect of increasing the width by a factor of ten. 

These results show that there is still good agreement between the full result 
and the spin correlation algorithm even when particles with fairly large widths are 
included, despite the potential problems. The result does not significantly depend on 
the choice of the spinor for the unstable particle and therefore we will use the second 
choice in order to avoid any dependence on the reference vector used to define the 
particle's spin. 



5. Tau Decays 

So far we have only discussed spin correlations in the SM for the top quark which 
is the only quark which decays before hadronization. However spin correlations can 
also be important for the decay of the tau. The tau will undergo a weak decay to a 
tau neutrino and either an electron/muon and its associated antineutrino or hadrons. 
Due to the relatively small mass of the tau only a small number of hadrons can be 
produced in its decay. 

The matrix elements for the leptonic decay of the tau and many of the hadronic 
decay modes are known, either theoretically or from experimental measurements of 
the distributions of the decay products. The helicity of the decaying tau can affect 
the distribution of the decay products. There is a package TAUOLA [32] which 
includes the matrix elements for the leptonic and many of the hadronic modes and 
includes the helicity of the decaying tau. 

In order to implement the algorithm we have suggested in full for tau decays 
we would need to completely rewrite the TAUOLA package to use the spin density 
matrices for the decays and return the decay matrix after performing the decays. 
Given the number of decay modes and the sophisticated treatment of the decays in 
TAUOLA this would be a major project and is beyond the scope of this study. 

However as TAUOLA makes use of the helicity of the decay tau to perform the 
decays we can use the diagonal entries of the spin density matrices, which give the 
probability the tau has a given helicity, to select the helicity of the decaying tau. 



6. Including Correlations with the Parton Shower 

Given the similarities between the algorithm we are using in order to include spin 
correlations from heavy particle decays and the algorithm of [5-8] for spin correlations 
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in the QCD parton shower it should be possible to produce an algorithm which 
includes both effects, i.e. correctly includes correlations from both heavy particle 
decays and in the QCD shower. 

In order to construct such an algorithm it is helpful to first review the algorithm 
of [5-8] for correlations in the QCD parton shower. The full algorithm is given for 
both forward and backward evolution in [7]. The algorithm proceeds in the following 
way: 

1. The momenta of the particles in the hard collision process are generated ac- 



cording to Eqn. 3.1. 



2. One of the outgoing partons is chosen at random and given a spin density 
matrix according to Eqn. |3.3| . 

3. The type of the branching, i ^ j + k, together with the azimuthal angle and 
momentum fraction z, are generated according to 

PlA:VS,A.(-,0)ViS;A^(->0)^U;^A,A^, (6.1) 

where Vi\ {z, 0) is the splitting function for the branching given the helicities, 
Aj, Xj, Afc of the partons, p\ _y is the spin density matrix of the parton entering 
the vertex and D-!^ ^ is the decay matrix for the partons leaving the vertex. 
As before in the first stage of the algorithm the decay matrix is taken to be 



4. One of the outgoing partons is selected to be developed and a spin density 
matrix calculated for this parton using 

"La; = t-WS:y,yS'' f^KK- (6-2) 

where the normalization 

iVpQCD = pLa;^4a.(^' <^)^A;t.A^(^' (6-3) 

is chosen so that the spin density matrix has unit trace. This procedure is 
repeated until the parton to be developed reaches the cut-off scale. As before 
the decay matrix for this parton is then set to ^AiA'- 

5. Another parton from the branching which produced the parton which reached 
the cut-off scale is now developed with a spin density matrix calculated for this 
parton using Eqn. |6.2| . When all the outgoing partons in a branching have been 
developed a decay matrix for that branching is calculated using 
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where the normahzation 

Ndqcb = V^X,^{z, <P)Vl^:,,,{z, <P)Di^yDl,, (6.5) 

is chosen so that the trace of the decay matrix is one. The other outgoing 
partons from the previous branching are now developed with a spin density 
matrix calculated using Eqn.|E]^ and the calculated decay matrix rather than 
the identity for those partons which have been developed. 

6. This procedure is repeated until the hard process is reached. Another outgoing 
parton is selected and a spin density matrix calculated for it using Eqn. ^]3| and 
the calculated decay matrices for those particles which have already been de- 
veloped rather than the identity. Step three of the algorithm is then performed 
for this particle. This procedure is repeated until all the outgoing particles in 
the hard process have been developed. 

This is the algorithm for the forward evolution of outgoing time-like partons. A 
similar algorithm [7] can be used to develop the incoming space-like partons. This 
algorithm starts the development with a decay matrix calculated according to 

Dk k' = —p'^^'M^n -x, X Ml,,,.y y TT V- (6-6) 

^ i Nj '^''^i ^I'^ji-^l---^"'' K-Kj;Xj^...X'„ J_J_ AfeAj, V / 

k=l,n 

The normalization 

k=l,n 

is again chosen such that the decay matrix has unit trace. We will not discuss the 
details of this backward evolution algorithm here, it is described in [7]. It is sufficient 
to note that when a space-like parton has been completely developed the algorithm 
returns a spin density matrix for the developed parton. So the full algorithm for spin 
correlations in the QCD shower is that the outgoing time-like partons are developed 
according to the algorithm described above and then using the decay matrices re- 
turned by this algorithm the incoming space-like partons are developed according to 
the algorithm of [7]. 

In order to treat the spin correlations the best possible approach would be as 
follows: 

1. The momenta of the particles in the hard collision are generated according to 
Eqn.|33. 

2. One of the outgoing particles is selected at random and a spin density matrix 
calculated for it according to Eqn.|3.3|. 
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3. If the particle is colourless and unstable step four is performed, if the particle 
is coloured step five is performed, otherwise the decay matrix for this particle 
is set to 5\.xi,, another particle selected and this step repeated. 

4. The decay mode of the particle is selected according to the branching ratios 
and the momenta of the particles produced in the n— body decay generated 
according to Eqn.|3.5|. 

(a) A particle from this decay is selected and a spin density matrix calculated 
for it using Eqn. |3.6| . 

(b) If this particle is unstable and colourless we perform step four. 

(c) If the particle is coloured step five is performed. 

(d) If the particle is stable and colourless its decay matrix is set to 5x-\i, 
and another particle from the decay which produced it is selected to be 
developed. Once all the particles in a given decay have been developed a 
decay matrix is calculated for the decay using Eqn. p.8| and another particle 
from the decay or branching which produced this particle selected to be 
developed. A spin density matrix is calculated for the development of 
this particle with the calculated decay matrix rather than the identity for 
those particles which have been developed. Step three is then performed 
for this particle. 

5. If the particle is coloured and has not reached the cut-off scale its branching, 
i — jTc, is generated according to Eqn. |6.1| . If the parton has reached the cut-off 
scale then: 

(a) If the particle is stable its decay matrix is set to 5\^y^ and another particle 
from the branching or decay which produced this parton selected to be 
developed. If all the particles in a branching have been developed a decay 
matrix for the parton which branches is calculated using Eqn.|6]^ and 
another particle from the decay or branching which produced it selected to 
be developed. Similarly if all the particles in a decay have been developed 
a decay matrix is calculated for the decay using Eqn. ^.Sl and another 
particle from the decay or branching which produced it selected to be 
developed. A spin density matrix is calculated for the development of 
this particle with the calculated decay matrix rather than the identity for 
those particles which have been developed. Step three is then performed 
for this particle. 

(b) If the particle is unstable step four is performed for the particle. 
If a branching occurs the following procedure is performed: 
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(a) One of the outgoing particles is selected to be developed and a spin density 
matrix calculated using Eqn. |6.2| . 

(b) Step five is then performed for this particle. 

6. Eventually this algorithm will return a decay matrix for one of the particles 
produced in the hard process. If there are still outgoing particles which have 
not been developed then another outgoing particle is selected to be developed 
and step two performed for this particle. 

7. If all the outgoing particles have been developed then one of the incoming par- 
ticles is selected and a decay matrix calculated using Eqn.|6.6|. The algorithm 
described in [7] is used to develop the incoming particle. This is repeated until 
all the incoming particles have been developed. 

As with both the QCD and heavy particle decay algorithms we have already 
described the decay matrix for any particle which has not been developed is taken 
to be 5x,x'. 

While this algorithm is the most efficient way in which to include the spin cor- 
relations it presents problems in generating the parton shower. This is because the 
parton shower algorithm needs to produce the QCD radiation from all the particles 
produced in the same hard collision or decay at the same time. There are several 
reasons for this: 

1. The scale used for the strong coupling constant in the parton shower algorithm 
depends on the energy of the collision, or mass of the decaying particle. 

2. In the angular-ordered parton shower used in HERWIG [33-35] the maximum 
angle for the emission of QCD radiation from a parton is determined by the 
direction of its colour partner, which is defined by the colour fiow in the hard 
process or decay. 

3. The parton- shower algorithm does not conserve energy and momentum. There- 
fore after the radiation from all the partons in a process has been generated 
the momenta of the resulting jets must be adjusted in order to ensure energy 
and momentum conservation. 

We need to consider an algorithm in which this is possible. While this algorithm 
will be less efficient in the way it implements the spin correlations it should prove 
much easier to interface with the QCD parton shower. The algorithm will proceed 
in the following way: 

1. First the algorithm for spin correlations in the both the initial- and final-state 
QCD parton shower is run as before. 
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2. Then if there are any unstable particles produced one of them is selected at 
random and a spin density matrix calculated in the following way: 

(a) The decay matrices for the other undecayed particles are set to ^a^a'- 

(b) We then use the matrix elements of all the vertices to produce the decay 
matrices of the particles in the hard process, apart from the particle which 
produced the one we have selected, and the spin density matrices for the 
incoming partons. This should be straightforward as these decay matrices 
will have already been calculated by the QCD spin correlation algorithm. 

(c) Using these decay matrices and Eqn. |3.3| we construct the spin density ma- 
trix for the parton which initiated the shower which produced the selected 
particle. 

(d) We then move down the branch that produced the particle we have se- 
lected using Eqn. |0| and the decay matrices which have already been 
calculated for the side-branches. 

(e) Eventually we will reach the particle we have selected with a spin density 
matrix calculated according to Eqn.|6.2|. 

3. The decay mode of this particle is selected according to the branching ratios 
and the momenta of the decay products generated according to Eqn. p75| . 

4. If any coloured particles are produced the QCD parton-shower algorithm is 
run. If the decay produces unstable particles we perform steps two and three 
for this decay, using Eqn.|3]^ rather than Eqn. ^TSj 

5. Eventually we will perform a decay in which only stable particles are produced. 
When this happens we set the decay matrices for the stable particles to ^a^a;. 
and use the decay matrices calculated by the QCD algorithm and Eqn.|3]^ to 
calculate a decay matrix for this decay. 

6. We then select another unstable particle produced in the previous decay, or 
the parton shower initiated by that decay, and repeat the procedure with the 
decay matrix calculated from Eqn. |3.8| for the particles for which the decays 
have been performed. If there are no further unstable particles to decay we use 
Eqn. |3.8| to calculate the decay matrix and move up the chain. 

7. Eventually we will reach a particle which was produced in the parton shower 
after the hard collision. We then select another unstable particle and repeat 
the process using the calculated decay matrix rather than the identity for those 
particles which have been decayed. 

8. This procedure is repeated until all the unstable particles produced in the 
parton shower initiated by the hard process have been decayed. 
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This procedure is more complicated as additional spin density and decay matrices 
have to be stored and recalculated but it should prove easier to implement the parton 
shower with this procedure. 

The implementation of either of the algorithms we have described in this sec- 
tion would involve a major rewrite of the HERWIG, or any other, event genera- 
tor. As there are currently projects under way to rewrite the current generation of 
FORTRAN event generators in C-|— I- [36,37] it would be pointless to undertake this 
rewrite. Hopefully, one of the algorithms we have described here can be implemented 
in the next generation of C-|— I- event generators. 

7. Conclusions 

We have shown that it is possible to construct an algorithm for spin correlations 
in a Monte Carlo event generator that both allows us to generate the production 
and decay of the heavy particles as separate steps and has a complexity which only 
grows linearly with the number of final-state particles. An additional advantage of 
this algorithm is that we only have to calculate the hehcity amplitudes for 2 — > 2 
collision processes and, at least for the processes we have studied, at most three body 
decays rather than the calculation of n— body matrix elements. 

This algorithm gives results that are in good agreement with the full calcula- 
tions of the n— body matrix elements for all the processes we have studied. In the 
appendix we have given the helicity amplitudes necessary to implement this algo- 
rithm in a Monte Carlo event generator. These results have been incorporated into 
the HERWIG event generator and will be available in the next version. 

Finally, we have proposed an algorithm to include both the spin correlations in 
heavy particle decays which we have considered here as well the correlations in and 
between jets which are already included in the HERWIG event generator. This algo- 
rithm, which would require a major rewrite of the current generation of FORTRAN 
event generators, can hopefully be included in the next generation of C-|— |- event 
generators which are currently being written [36, 37] . 
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A. Spinor Conventions 

In general we will use the conventions of [38, 39] for the spinors in the helicity am- 
plitudes for the matrix elements. However it is more convenient for us to define 

S\{Pl,P2) =Ux{pi)u_x{P2), (A.l) 

rather than the s and t functions of [38] . 

We will use the following notation for massive fermions, 

u{p,X)= -j^={]^+m)u_x{l), (A.2a) 

v{p, A) = -^yi^ (/-H «-a(0, (A.2b) 

where p is the four-momentum of the fermion, m is its mass, A is the spin of the 
fermion and £^ is a light-like four- vector. This choice for the massive fermions both 
reproduces the standard spin sums for Dirac particles 

y^^uxip)ux{p) = (/+m) , (A.3a) 

A 

J2'"x{p)vx{p) = (y-m), (A.3b) 

A 

and the additional spin sums required in SUSY models which have Majorana fermions 
[40] 

A 

J2 ul{p)vx{p) = i^-m) , (A.4b) 

A 

J2 vl{p)ux{p) = i^+m) , (A.4c) 

A 

Yl MpWxiP) = i^-m) C^. (A.4d) 

A 

Given this choice it is useful to define 

where — Q and £^ = 0. In this notation the spin vector for the particle is given by 

m 

The choice of the spin vector is completely arbitrary in our case provided that we are 
consistent, i.e. we use the same choice for the production and decay of the particle. 
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In the expressions for the amphtudes it will often be convenient to express the 
results in terms of the function 



F{\i,pi,k,m,\2,P2) = ux^ipi) (f:+m) ux^{p2), (A. 7) 

where pi^2 are massless four vectors, Ai^2 are helicities and k is an arbitrary four- 
vector. This function can be evaluated in terms of spinors as follows [39] 

F(A, pi, k, m, -A, P2) = msx{pi,P2) (A. 8a) 

F{X,pi,k,m, X,p2) = sx{pi,k)s^x(k,p2), (A. 8b) 

where 

We will use the same notation as [38] for the polarization vectors of the massless 
gauge bosons 

e^(p) = -±=ux{p)Yuxil), (A.IO) 

where p is the four-momentum of the gauge boson and / is an arbitrary light-like four- 
vector which is not coUinear to p. The choice of this light-like vector corresponds to 
the making a choice of gauge and therefore we can either chose l'^ in order to simplify 
the calculation or vary it in order to test the gauge invariance of the results. 

The choice of the representation for the polarizations of the massive gauge bosons 
is discussed in [38]. However, as this representation basically replaces the gauge boson 
with two massless fermions, into which the boson decays, we will not include any 
on-shell massive gauge bosons but replace them with their decay products. 



B. Couplings 

In this section we will give the couplings for both the SM and MSSM in a format 
suitable for use in the helicity amphtudes for both the production and decay of the 
particles. 

In order to specify the matrix elements we need the couplings for the vertices 
involved. In general we have taken the form of a vector-fermion-fermion vertex to be 

iaV^A- (B.l) 

The couplings of the gauge bosons to the Standard Model fermions are given in 
Table|^ where e is the magnitude of the electron's charge, g the weak coupling and 
6w the weak mixing angle. 
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C^Bff 



tWff 



Bff 



Bff 



'^Zff 



cos 9w 



cos 6 



(t^ - Cf sin^ ^ 



Table 3: Couplings of the SM fermions to the gauge bosons, where is the weak isospin 
of the fermion and ef is its electric charge in units of the magnitude of the electron's charge. 
The isospin partner of the fermion f is denoted by f. Vf[i is the relevant element of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix. 





'^Bxx 


«Bxx 




{U:,U,, + \U*,U,2 - S., sin^ 0w) 




«ZxOx? 






"Wx9x- 







Table 4: Couplings of the gauge bosons to the electroweak gauginos. The mixing matrices 
for the neutralinos, N, and charginos, U and V, are defined in the notation of [40-42]. 



In the MSSM there are a lot more couplings which must be specified and this 
is complicated by the mixing of the particles involved. The physical electroweak 
gauginos are mixings of the interaction eigenstates. There are mixing matrices: 
for the neutralinos; and U and V for the charginos. In general our conventions for the 
mixing of the gauginos follow those of [40-42]. In performing the diagonalization to 
go from the interaction to mass eigenstates there is a choice: we can either take the 
mixing matrices to be complex, in which case the masses are positive; or the mixing 
matrices to be real and let the masses be positive or negative, in which case we must 
redefine the field x ~^ 75X if mass is negative [40-42]. As the calculations will 
be implemented numerically it is more convenient to deal with real mixing matrices 
and we must therefore keep track of the sign of the gaugino masses and change the 
sign of certain couplings. We will denote the sign of the mass of the ith neutralino 
by ej and the zth chargino by 77,. The couplings of the gauge bosons to the gauginos 
of the MSSM are given in Table^. 

In general there can be mixing between all the sfermions with the same electric 
and colour charges. However, we will only consider mixing between the left and right 
scalar partners of the same fermion. We use the same convention as the HERWIG 
event generator for this left/right mixing, this is discussed in [16,26,27]. In this 
notation the mixing matrix for the squarks is Q^p, where i = 1 — 6 for the d, u, s, c, 
b and t squarks, a is the left/right eigenstate and (3 is the mass eigenstate. Similarly 
the mixing matrix for the sleptons is L^^, where i = 1 — Q for the e, u^, /i, z/^, r, 1/^ 
sleptons, a is the left/right eigenstate and /3 is the mass eigenstate. 
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'^Bff* 




'^Bff* 


















V2^1" 




—iecf 







Table 5: Couplings of the sfermions to the electroweak gauge bosons. These couplings 
are taken from [40]. The flavour of the sfermion is denoted by i and a,f3 gives the mass 
eigenstate. 



obhh 




Qbhh 




aw±H±ho 


=F| cos(/9 — a) 


aw±H±AO 


*2 


aw±H±HO 


±f sin(/3 - a) 




j9cos{/3-a) 
2 cos 9w 




^gsin{l3-a) 
2 cos 6w 







Table 6: Couplings of the Higgs bosons to the electroweak gauge bosons. These couplings 
are taken from [43,44]. The angle a is the mixing between interaction and mass eigenstates 
of the two neutral scalar Higgs bosons. 



The form of the boson-scalar-scalar vertex is taken to be 

ia{p + pY, (B.2) 

where the momentum p is taken to be flowing into the vertex and p' is taken to be 
flowing out of the vertex. The couplings of the sfermions to the electroweak gauge 
bosons are given in Table|^ and of the Higgs bosons to the gauge bosons in Tablej^. 
The form of the scalar-fermion-fermion vertex is 

ia^Px. (B.3) 

The couplings of the neutralinos to the sfermions are given in Table 0, of the charginos 
to the sfermions in Table ^ and of the gluino to the squarks in Table ^. In order to 
define the couplings of the sfermions to the neutralinos it is useful to define the 
following functions 

5jo = -Nieef + N'^a+^-^, (B.4a) 
for the couplings of the gaugino part of the neutralinos to the fermions. The functions 



for up — type quarks and neutrinos, 



2Mw sin /3 

^fx? = { (B-5) 



2Mw cos/3 down — type quarks and charged leptons. 
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) 




) 








-y2(L*2„i/*-o + I^l„5-.o) 





Table 7: Couplings of the neutralinos to the sfermions and fermions. The flavour of the 
sfermion is denoted by i and a gives the mass eigenstate. 



give the couplings of the fermions to the Higgsino part of the neutralinos. 

The couplings of the Higgs boson of the MSSM to the electroweak gauginos are 
given in Table ^ where 



Q^ 



Q 



iR 



Qll 



Vi 
V2 



j2, 



v2 



1 



N^3UJl - -^Uj2 {Ni2 + AT.i tan^w-) 



Q'L = i [N^3 {N,2 - N^i tan Ow) + iV,3 (iV.2 - Na tan Ow)] 



[Nii {Nj2 - Nji tan Ow) + Nj^ {Ni2 - Na tan Ow)] 



(B.6a) 
(B.6b) 

(B.6c) 

(B.6d) 

(B.6e) 
(B.6f) 



taken from [43]. 

The couplings of the Higgs bosons to the Standard Model fermions are given in 
Table0. 

We also need the couplings of the MSSM Higgs bosons to gauge boson pairs. 
The coupling has the form iag'^'^ with the couplings gM-\j^ sin(/5 — a) for the lightest 
Higgs boson to a pair of either W or Z bosons. The coupling of the heavier scalar 
Higgs boson is (y'Mw cos(/3 — a). 



C. Production Matrix Elements 



There are a large number of both SM and MSSM processes implemented in most 
general purpose Monte Carlo event generators. In the case of the Standard Model 
most of these processes do not involve the production of the top quark or tau lepton 
which are the only fermions for which the spin correlations are relevant. There are 
a number of other processes involving three or four particles in the final state, for 
example production of a gauge or Higgs boson in association with a tt pair which 
have much smaller cross sections. 



37 









%+u d 


V2M^cos/3"^diC//2Qi„ 












) 













^ ( JT* T2i-l 1 rr* r2j-l\ 



Table 8: Couplings of the sfermions to the charginos. The generation of the sfermion is 
given by i and a is its mass eigenstate. 



flgqq 


gqq 






9sV2Qf^ 


-9sV2Qf^ 




9sV2Qf-' 


-9sV2Qf-' 



Table 9: Couphngs of the squarks to the gluino. The flavour of the squark is given by i 
and a is the mass eigenstate of the squark. 









«hox?X° 


g [Q''j sin a + S'/j cos a) 


g {Q'lj* sin a + S'^* cos a) 


«Hox?X? 


g i^S'lj sin a — Q"j cos a) 


g [S'j* sin a — Q'-* cos a) 


«Aox?X? 


-ig {Q'r sin f3-S^j cos 


ig {Q'^t ^in 13 -S';t cos 13) 


%x+x7 


g {Qij sin a — Sij cos a) 


g {Q*j^ sin a — S*^ cos a) 


^Hox+x7 


—g {Qij cos a + Sij sin a) 


—g {Q*i cos a + 5**^ sin a) 


^Aox+X" 


-ig (Qij sin (3 + Sij cos/3) 


ig {Q*jiSin 13 + S*jiCos 13) 


^H-x?x7 


-gQ'^ sin p 


-gQfj cos (3 



Table 10: Couplings of the MSSM Higgs bosons to the gauginos. 





a+- 

"Hff 


'^Hff 


OHff 


"Hff 


'^Hff 




gniui cos a 






pmuj sin a 




2Mw sin f3 


2Mwsin/3 


"hOdidi 


gm^^ sm a 


"h'M.d, 


'^H"did, 


3m;j^ cos a 


^HOd.d, 


2Mw cos 13 


2Mw cos /3 




gmi^ sin a 
2Mw cos /3 






gm^^ cos a 
2Mw cos /3 


a+ - 




. gniu - cot /3 
^ 2Mw 




'^H+Uid, 


flmd, tan /3 
V2Mw 


gmu^ cot P 
\/2Mw 




. gm^ . tan /3 
* 2Mw 


"AOdidi 




gm^^ tan /3 
\/2Mw 





".V'lj:, 


. f/w^ . tan 
' 2Aiw 











Table 11: Couplings of the MSSM Higgs bosons to the Standard Model fermions. 
A large number of the MSSM production processes involve the production of a 

pair of scalar particles, and therefore we do not need the matrix elements for these 
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processes.^ This leaves a reasonably small number of 2 ^ 2 processes for which we 
must calculate the matrix elements. 

In this section we will first present the helicity amplitude expressions for all the 
Feynman diagrams which occur in the 2^2 processes we are considering. These 
amplitudes can then be combined to give all the matrix elements we will need. In all 
the matrix elements Pi and rrii are the four-momentum and mass of the ith particle, 
respectively. The vector li is used to define the direction of the particle's spin and 
Aj is the helicity of the ith particle. 

C.l Feynman Diagrams 
C.1.1 Diagram 1 




Figure 18: Feynman diagram for f(pi)f(p2) f'(P3)f'(P4) via s-channel gauge boson 
exchange, where f and f can be any fermion. 

The Feynman diagram for ff f'f via s-channel gauge boson exchange is shown 
in Fig.|IH[ The form of the first and second vertices are taken to be ia^'y^Px and 
ib^'y'^Px. The amplitude is given by 



M = - , . : (c.i) 



VP3 ■ Im ■ h S-MI + iVbMb 

+b^^ F(-A3,/3,p3,m3, -Xi,p2)F{-\i,pi,p4, -m4, -A4, /4 

where s = {pi +P2Y, Mb is the mass of the exchanged boson and F^ is the width of 
the exchanged boson. The masses of the incoming particles have been neglected. If 
we wish to consider both the outgoing particles to be fermions rather than a fermion 
and an antifermion, in for example gaugino pair production, the sign of A4 should be 
changed. 

C.l. 2 Diagram 2 



The Feynman diagram for ff ^ f'f via t-channel scalar exchange is shown in Fig. |T9 
We will take the form of the upper vertex to be ia^P\ and the lower vertex to be 



^The matrix elements used for these processes in HERWIG are given in [16,28]. 



39 



1 > ^ > 3 



I 



6^ 



Figure 19: Feynman diagram for f(pi)f(p2) f'(P3)f'(P4) via t-channel scalar exchange, 
where f and f can be any fermion. 

ih^P\. The matrix element is given by 

„AiL-A2 1 

M = ^ —— (C.2) 

2 VP3 ■ hPA -ht-Ml 

F(-A3, /3,p3, ma, Ai,pi)F(A2,P2,P4, -m^, -XaJa), 

where i = [pi — and M$ is the mass of the exchanged scalar. The incoming 
particles have been taken to be massless. As before if we wish to regard particle four 
as a fermion rather than an antifermion the sign of its helicity should be changed. 

C.1.3 Diagram 3 

1 ^ ^ ^ 4 



Figure 20: Feynman diagram for f(pi)f(p2) f'(P3)f'(P4) via u-channel scalar exchange, 
where f and f can be any fermion. 

The Feynman diagram for ff f'f via w-channel scalar exchange is shown in 
Fig.^. As for the previous diagram we will take the form of the upper vertex to be 
ia^Px and the lower vertex to be ih^P\. The amplitude is given by 



F(A4, k,p4,m^, \i,pi)F{p2, X2,P3, -ms, A3, 13), 



where ii = {pi — ^4)^, M$ is the mass of the exchanged sfermion and we have assumed 
the incoming particles are massless. The sign of A4 should be changed if we wish to 
consider both outgoing particles as fermions. 
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Figure 21: Feynman diagram for g(pi)g(p2) f(P3)f(P4) via t-channel fermion exciiange. 



C.1.4 Diagram 4 



The Feynman diagram for gg — * ff via t-channel fermion exchange is shown in Fig. 21 



In this case the form of the gluon-quark-quark vertex is taken to be —igs'^^-, where Qs 
is the strong couphng. The colour matrices are not included as it is easier to handle 
the colour sums/averages for the diagrams separately. This would also allow us to 
use the same matrix elements for fermion- ant if ermion production in photon-photon 
collisions by replacing the strong coupling with the electric charge. The matrix 
element is given by 



M 



9l 



1 



(C.4) 



2 VPi -Im- kP3 ■ hPi ■ h t 
[ F(-A3,/3,p3,m3, Ai,/i) 

{ F{ Ai,Pi,P2 -P4,"^3, ^2,h)F{ \2,P2,P4,-mi,-\^,k) 

+F{ \l,Pl,P2 -P4,ms, -X2,P2)F{-X2j2,P4, -"^4, -A4,/4)} 
+F{-\3,l3,p3,m3, -Xi,pi) 

{ F{-Xi,li,p2 - P4,m3, X2,h)F{ X2,P2,P4,-m4,-X4,,U) 
+F{-Xi,h,P2 -P4,m3, -X2,P2)F{-X2,l2,P4., -"^4, -A4,/4)}] • 

As before if we wish to regard the outgoing antifermion as a fermion, in for example 
gluino pair production, then we must change the sign of A4. 

C.1.5 Diagram 5 




Figure 22: Feynman diagram for g(pi)g(p2) f(P3)f (P4) via w-channel fermion exchange. 



The Feynman diagram for gg ff via w-channel fermion exchange is shown in 
Fig.E^. As with the previous diagram we have not included the colour matrices in the 
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expression for the helicity amplitude. In this case the form of the gluon- quark- quark 
vertex is taken to be —igs'^^- The matrix element is given by 



M 



9l L_ 

2a/Pi ■ hP2 ■ hp?, ■ hPA ■ h u - ml 



(C.5) 



[ F{-\3,k,ps,ms, \2,l2) 

{ F{ X2,P2,Pi - Pi,ms, Xi,h)F{ Ai,pi,p4, -m4, -A4,/4) 
+F{ X2,P2,Pi -P4,rns, -Ai,pi)F(-Ai, Zi,p4, -m4, -A4,Z4)} 

+F(-A3,/3,P3,m3, -A2,P2) 

{ F{-X2,l2,Pi - P4,m3, Xi,h)F{ Xi,pi,p4,-m4,-X4,,U) 
+F{-X2,l2,Pi -P4,m3, -Ai,pi)F(-Ai,/i,p4, -m^, -A4,/4)}] • 



As before if we wish to regard the outgoing antifermion as a fermion, in for example 
gluino pair production, then we must change the sign of A4. The sign of this matrix 
element has been changed in order to simplify the combination of the matrix elements 
with the colour factors. 



C.1.6 Diagram 6 




Figure 23: Feynman diagram for g(pi)g(p2) f(P3)f(P4) via s-channel gluon exchange. 



The Feynman diagram for gg — > ff via s-channel gluon exchange is shown in 
Fig.^. In this case the form gluon-quark-quark vertex is taken to be —igsj^ and 
the triple gluon vertex to be igs [{pi —piVg"^^ + {pi + '^P2)"^g^^ — (2pi + P2Y g°''^'^^, 
where P2il3 and pi + P2,'J are the momenta and Lorentz indices of the three 

gluons. The colour matrices are not included in the expression for the amplitude 
in order to allow us to use the same matrix element for both quark and gluino 



42 



production. The matrix element is given by^ 

1 

M = (C.6) 
4vPi ■ hP2 ■ I2P3 ■ hP4 ■ h s 

[( 5ai-A2SAi(Pi,^2)S-Ai(P2,^i) + 5aiA2-SAi(Pi,P2)S-Ai(^2,^i)) 

( F(-A3,/3,p3,m3, A,pi)F(A,pi,p4, -m^, -K,h) 

-F{-Xs, I3, J93, m3, A, P2)F{\ P2, P4, -"^4, -A4, U)) 
+2SAi(Pi,P2)S-Ai(P2,^i) 

( F(-A3,/3,P3,m3, X2,h)F{ A2,P2,P4, -"^4, -A4, ^4) 
+F(-A3,/3,P3,m3, -A2,P2)i^(-A2,^2,P4, ""^4, -A4, ^4)) 
-2Sa2(P2,Pi)S-A2(Pi,^2) 

( F{-\3,l3,P3,ms, Ai,/i)F( Ai,Pi,P4, -"^4, -A4,^4) 

+F(-A3,/3,p3,m3, -Ai,pi)F(-Ai,/i,p4, -m4, -A4,/4))] . 

As before if we wish to regard the outgoing antifermion as a fermion, in for example 
gluino pair production, then we must change the sign of A4. 

The expressions we have given for gg ff via t- and w-channel fermion exchange 
and s-channel gluon exchange are given in a general gauge. While we used these 
results to check the gauge invariance of the various production processes, in HERWIG 
we made the gauge choice li = p2 and I2 = Pi which simplifies the expression given in 
Eqn. |C6| for the s-channel gluon exchange which involves the triple gluon vertex [38]. 

C.1.7 Diagram 7 



t 

I 
I 



Figure 24: Feynman diagram for f(pi)g(p2) x(P3)f(P4) via t-channel sfermion ex- 
change. 

The Feynman diagram for fg — > xf via t-channel sfermion exchange, where x 



can be any of the gauginos, is given in Fig.^. The form of the fermion-gaugino- 
scalar coupling is ia^P\ and the form of the scalar-scalar-gauge boson vertex is 
—Ws{p + p')^- As before the colour matrices are not included in the amplitude. 
The helicity amplitude is given by 

9s Q-^^ 

M = ' =j ji^(-A3,/3,P3,"^3, Al,Pl)F(A2,P2,P4,0, A2,/2), (C.7) 

V 2P2 ■ ^2P3 -ht-ml 



*The polarization A should be summed over. 
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where the incoming fermion is assumed to be massless. 



C.1.8 Diagram 8 




Figure 25: Feynman diagram for f(pi)g(p2) x(P3)f(P4) via s-channel quark exchange. 



The Feynman diagram for fg yi via s-channel quark exchange is shown in 
Fig.^. The form of the fermion- fermion-gauge boson vertex is —igsl^ and the form 
of the fermion-fermion-scalar vertex is ia^P\. Again x can be any of the gauginos 
and the colour matrices are not included in the amplitude. The helicity amplitude 
is given by 



M 



9s 



1 



\/2p2 ■ hps ■ h S 

[ Sx^x^a^^sx^ {pi,p2)F{-X3, /g, p3, ms, A2, Pi)s_a2 (h, Pi) 

+5Ai-A2a"^'-5A2(p2,Pl) ( F(-A3,Z3,P3,m3, -A2,Pl)s_A2(Pl,^2) 

+ F(-A3,/3,P3,m3, -A2,P2)S_A2 iP2,l2))] 



(C.8) 



where the incoming fermion is assumed to be massless. 



C.1.9 Diagram 9 



Figure 26: Feynman diagram for f(pi)g(p2) x(P3)f (P4) via u-channel gaugino exchange. 

The Feynman diagram for fg — > xf via w-channel gaugino exchange is shown in 
Fig.^, due to the couphngs this diagram only occurs when the outgoing gaugino 
is the gluino. The form of the fermion-fermion-scalar vertex is taken to be ia^Px 
and the form of the gaugino-gaugino-gauge boson vertex is taken to be igsl^- As 
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before the colour matrices are not included in the amplitude. The helicity amplitude 
is given by 



M = - , , , (C.9) 
V 2p2 ■ hps -hu-mi, 

[ F{-\3,l3,p3,m3, X2,h)F{ X2,P2,Pl - Pi,m3,\i,pi) 
+F{-\3,l3,p3,m3, -\2,P2)F{-\2,l2,Pl " P4, "^3, Ai, Pi)] , 

where the incoming fermion is assumed to be massless. 

As with the diagrams for gg —>■ ff the amplitudes were implemented with a 
general choice of the gauge vector for the incoming gluon and checked for gauge 
invariance. However it is again convenient to make the gauge choice I2 = Pi in order 
to simplify the Feynman diagrams. 

C.1.10 Diagram 10 

1 " ^ ► 3 



I 
I 



Figure 27: Feynman diagram for f(pi)g(p2) x(P3)f*(P4) via t-channel sfermion ex- 
change. 

The Feynman diagram for fg —>■ xf * via t-channel sfermion exchange is shown in 
Fig.^. The outgoing gaugino can be any of the gauginos. The form of the fermion- 
fermion-scalar vertex is taken to be ia^Px and the scalar-scalar-gauge boson vertex 
is taken to be —igs{p + p')^, where we have not included the colour matrices in the 
amplitude. The helicity amplitude is given by 

M = ' =j ^F{\i,pi,p3,-m3,\3,h)F{X2,P2,PA,0,X2,P2), (C.IO) 

V 2P2 ■ hps ■ 1st -mi 

where the incoming antifermion is assume to be massless. 
C.1.11 Diagram 11 

The Feynman diagram for fg xf* via s-channel fermion exchange is shown in 
Fig.^. The outgoing gaugino can be any of the gauginos. The form of the fermion- 
fermion-gauge boson vertex is taken to be —igs'y'^ and the form of the fermion- 
fermion-scalar vertex is taken to be ia^P\. As before the colour matrices are not 
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Figure 28: Feynman diagram for f(pi)g(p2) x(P3)f*(P4) via s-channel quark exchange. 

included in the amplitude. The helicity amplitude is given by 

M = , - (C.ll) 

+ S_A2 (^2, P2)F{\2, P2,P3, -nis, A3, I3)) 

+5x^-\2a^^sx^ (pi, l2)sx2 {p2,Pi)F{-X2, pi, Ps, -ms, A3, h 
where the incoming antifermion is assumed to be massless. 
C.1.12 Diagram 12 



Figure 29: Feynman diagram for f(pi)g(p2) — > x(P3)f*(P4) via u-channel gaugino ex- 
change. 

The Feynman diagram for fg xf* via w-channel gaugino exchange is shown in 
Fig.^, due to the couplings in this case the outgoing gaugino must be the gluino. 
The form of the fermion-fermion-scalar vertex is taken to be ia^Px and the form of 
the gaugino-gaugino-gauge boson vertex is taken to be igs'j'^- The colour matrices 
are not included in the expression for the amplitude. The helicity amplitude is given 

by 

M = - , , (C.12) 

[ F{\i,pi,pi-p4,-m3, X2,l2)F{ \2,P2,P3,-fn3,\3,l3) 
+F(Ai,Pi,Pi -P4, -1713, -A2,P2)^(-A2,/2,P3, ""^3, A3, /3)] , 

where the incoming antifermion is assumed to be massless. 
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As with the diagrams for gg —>■ ff and fg — > xf the amphtudes were implemented 
with a general choice of the gauge vector for the incoming gluon and checked for 
gauge invariance. However it is again convenient to make the gauge choice I2 = pi 
in order to simplify the Feynman diagrams. 



C.1.13 Diagram 13 



Figure 30: Feynman diagram for f(pi)f'(p2) — > f(p3)f'(p4) via i-channel gauge boson ex- 
change. 

The Feynman diagram for ff — > ff via t-channel gauge boson exchange is shown 
in Fig. This diagram occurs in single top quark production via t-channel W 
exchange. The helicity amplitude is 

M = ^ l^^^^ ^L = [ (5ai-A2^(-A3,^3,P3,"^3, A2,P2)sA2(P4,Pi) (C.13) 

t - Vp3 ■ ^3 

+'^AiA2^(-A3, ^3, P3, "^3, -A2, P4)S-A2 iP2, Pi)] , 

where in addition to neglecting the masses of the incoming particles the mass of 
particle four, which will be a light quark in single top production, has also been 
neglected. 



C.1.14 Diagram 14 



Figure 31: Feynman diagram for f(pi)f (P2) — ^ f(P3)f (P4) via i-channel gauge boson ex- 
change. 

The Feynman diagram for ff — ff via t-channel gauge boson exchange is shown 
in Fig. This diagram occurs in single top quark production via t-channel W 
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exchange. The hehcity amphtude is 
t - Vp3 ■ h 

+5aiA2^(-A3, h,P3, ms, -A2, P2)S-A2 (P4, Pi)] , 

where in addition to neglecting the masses of the incoming particles the mass of 
particle four, which will be a light quark in single top production, has also been 
neglected. 

C.1.15 Diagram 15 



1 



2 




3 



4 



Figure 32: Feynman diagram for f(pi)f'(p2) f(P3)f'(P4) via t-channel gauge boson ex- 
change. 

The Feynman diagram for ff if via t-channel gauge boson exchange is shown 
in Fig. This process occurs in single top quark production via t-channel W ex- 
change. The helicity amplitude is 

M = ^i^'^" Sx,-X2Sx,{pi,P2)F{ X2,Pi,P3,-m3,-\3,k) {C.15) 

t - VP3 ■ h 

+^AiA2SAibl,P4)i^(-A2,P2,P3, -"^3, -A3, ^3)] , 

where in addition to neglecting the masses of the incoming particles the mass of 
particle four, which will be a light quark in single top production, has also been 
neglected. 

C.1.16 Diagram 16 

The Feynman diagram for ff —>■ ff via t-channel gauge boson exchange is shown 
in Fig. This diagram occurs in single top quark production via t-channel W 
exchange. The helicity amplitude is 

= X .^2 ' " / ^ [ ^X^-X2SXi{Pl,Pi)F{ A2,P2,P3,-"^3,-A3,^3) (C.16) 

t - VPs ■ h 

+5AiA2-SAi(Pl,P2)i^(-A2,P4,P3, -"^3, -A3, ^3)] , 

where in addition to neglecting the masses of the incoming particles the mass of 
particle four, which will be a light quark in single top production, has also been 
neglected. 
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Figure 33: Feynman diagram for f(pi)f'(p2) f(P3)f'(P4) via t-channel gauge boson ex- 
change. 

C.2 Matrix Elements 

We can now combine the amphtudes given in the previous section to give the matrix 
elements for the various processes. 

The matrix elements for the 2^2 Standard Model top quark and tau lepton 
production processes are given in Table p!^. It is easiest to extract the colour matrices 
from the diagrams and perform the traces separately. This leads to matrix elements 
for the different colour flows in the process and colour factors of the squares of 
the helicity amplitudes for the individual colour flows and the interferences between 
them. The colour flows for gg — >■ tt are discussed in Section W^. 



The matrix elements for the 2^2 MSSM electroweak gaugino pair, elec- 
troweak gaugino and gluino, and gluino pair production processes are given in Ta- 



bles 0, |1J and |I5|, respectively. The discussion in Section on the colour flows also 
applies to qq gg which has the same colour structure as gg — > tt, after crossing. 

The situation is slightly different for gg gg. As with qq gg the s-channel 
diagram which again involves the triple gluon vertex contributes to both colour flows. 
In this case the colour matrices which occur in the s-channel diagram can be rewritten 
using the Jacobi identity 

j^abe j^ecd yfece jaed jbed jace j^-^^ 

where a is the colour of the first incoming gluon, h is the colour of the second incoming 
gluon, c is the colour of the first outgoing gluino and d is the colour of the second 
outgoing gluino. This allows us to write the colour matrices for the s-channel diagram 
in terms of those for the t- and w-channel diagrams which means we can write the 
s-channel diagram as a piece which contributes to the t-channel colour flow and one 
which contributes to the w-channel colour flow. 

The matrix elements for the sfermion-gaugino production processes are given in 



Table 16. The colour flows of the gluino-squark production processes are the same. 



after crossing, as those in gg — > tt. 

D. Decay Matrix Elements 

In the Standard Model, for the processes we are considering, there are only two 
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Process 


NcF 


Colour Factors 


ID 


CF 


VP 










Cii 
















ff nn 


1 


1 






1 
1 


1 
1 


7 
Z 


— ecf 

ZjII 


-ecn 


qq' ^ f f ' 


1 


1 






1 


1 






- 

VV II 


qq tt 


1 


liVf 






1 


1 


O' 


1 


1 


gg^tt 


2 


1 

4Nc 


1 


1 

4Nc(N^-l) 


4 


1 


t 
















5 


2 


t 
















6 


1 


g 
















6 


2 


g 






bu — s> td 


1 


1 






13 


1 


w 


- 

VV t D 


'^Wdu 


bd tu 


1 


1 






14 


1 


w 


"wtb 


'^Wdu 


bd tu 


1 


1 






15 


1 


w 


"wtb 


'^Wdu 


bu td 


1 


1 






16 


1 


w 


"wtb 


^Wdu 



Table 12: Matrix elements for the Standard Model production processes. The number of 
different colour flows is given by Nqf, ID gives the type of diagram, CF gives the colour 
flow for a given diagram and VP gives the virtual particle exchanged in the s-, t- or u- 
channels depending on the diagram. The isospin partner of a fermion is denoted with a 
prime and n can be any fermion. The colour factor given for ff — ^ nh is for tau production 
in lepton collisions or electroweak production of top quark pairs in qq annihilation. The 
corresponding colour factor is 1 /Nc for tau production in qq annihilation. The colour factor 
given for ff production is for the production of tb or bt, the corresponding colour factor 
for t~Vt or t~^Vt production is l/N^. In the single top quark production processes the up 
and down quarks can be any of the quarks with the same weak isospin. 

decay modes, i.e. the decays of top and antitop quarks via a W boson. In the MSSM 
however, there are a large number of decay modes which must be calculated. As 
before we will only calculate those decay modes for which there are fermions and the 
use of the spin correlation algorithm becomes important. 

Despite the large number of decay modes there arc only a small number of 
Feynman diagrams. We have calculated these with arbitrary couplings which can 
then be specified for particular decay modes. We will first give the helicity amphtudes 
for these processes and then the diagrams and couplings involved for particular decay 
modes. 

D.l Two Body Decay Feynman Diagrams 

There are only five two body decay processes for which we need the helicity ampli- 



50 







Colour Factors 


ID 


CF 


VP 






Process 
























Cii 


















1 


1 


- 


- 


1 


1 


z 


- 

Zff 


^Zx^X? 












2 


1 




- 


"xOf.f 












3 


1 






*— A 

^*"x?Lf 


ff xtxj 


1 


1 






1 


1 


7 


— ecf 














1 


1 


Z 


- 

Zff 


Zxfx, 




only u-type quarks 


2 


1 


f' 

a 




*-A 




only d-type quarks/leptons 


3 


1 


r 

a 




*-A 




1 


1 

A^c 






1 


1 


w 


"Wud 


"wx^x" 












2 


1 


da 




X?d„d 












3 


1 


Ua 


■? x"UaU 


X, ""d 


du 


1 


1 






1 


1 


w 




"wx^x" 












3 


1 


Ua 


xr^^d 














2 


1 


da 




X, d^u 



Table 13: Matrix elements for the Minimal Supersymmetric Standard Model electroweak 
gaugino pair production processes. The number of different colour flows is given by Nqf, 
ID gives the type of diagram, CF gives the colour flow for a given diagram and VP gives 
the virtual particle exchanged in the s-, t- or li-channels depending on the diagram. The 
colour factors given for x'^x'^ &iid X^X~ production are for incoming leptons, the colour 
factor is 1/Vc for incoming quarks. The mass eigenstate of the sfermion is given by a and 
where the sfermions appear in u- and t-channel propagators the mass eigenstates should 
be summed over. The isospin partner of a fermion, or sfermion, is denoted by a prime. An 
identical particle symmetry factor of one half must also be included for the production of 

A.I • 



tudes and most of these processes are relatively simple as they involve scalars. The 
four-momentum, mass, helicity and reference vector used to define the ith particle's 
spin are pi, rrii, \i and k, respectively. 



D.1.1 Diagram 1 



The decay of a fermion to a fermion and a scalar is shown in Fig. 0. This process 
is common in SUSY models involving the decays of the gauginos to either a fermion 
and an antisfermion or to a gaugino and a Higgs boson. As before we will take the 
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NcF 


Colour Factors 


ID 


CF 


VP 






Process 




























C12 














1 




- 


- 


2 


1 


qa 




"gqaq 












3 


1 


(In 




X,qaq 




1 


(^1-1) 
2JV2 






2 


1 


da 




gdad 












3 


1 


Ua 


"gu„u 


X, Uad 




1 


(iv|-i) 

2JV2 






2 


1 


U„ 




gUaU 












3 


1 


da 


gdad 


xrdau 



Table 14: Matrix elements for the Minimal Supersymmetric Standard Model electroweak 
gaugino production in association with a gluino. The number of different colour flows is 
given by Nqf-, ID gives the type of diagram, CF gives the colour flow for a given diagram 
and VP gives the virtual particle exchanged in the s-, t- or u-channels depending on the 
diagram. The mass eigenstate of the squark is given by a and when the squarks appear 
in u- and t-channel propagators the mass eigenstates should be summed over. The isospin 
partner of a fermion, or sfermion, is denoted by a prime. 



Process 


NcF 


Colour Factors 


ID 


CF 


VP 






Cii 


C22 




qq gg 


2 


(iV2-l)2 

8JV^ 


(iV2-l)2 


mi 


1 


1 


g 


1 


1 


1 


2 


g 


-1 


1 


2 


1 


qa 


gqaq 


gq^q 


3 


1 


qa 


a^- 
gqaq 


gqc«q 


gg^ gg 


2 


2(iV2-l) 


2(JV^-1) 




4 


1 


g 






5 


2 


g 






6 


1 


g 






6 


2 


g 







Table 15: Matrix elements for the Minimal Supersymmetric Standard Model gluino pair 
production processes. The number of different colour flows is given by Ncf, ID gives 
the type of diagram, CF gives the colour flow for a given diagram and VP gives the 
virtual particle exchanged in the s-, t- or ti-channels depending on the diagram. The mass 
eigenstate of the squark is given by a and where the squarks appear in u- and f-channel 
propagators the mass eigenstates should be summed over. The isospin partner of a quark, 
or squark, is denoted by a prime. The identical particle symmetry factor for gg production 
has been included in the colour factor. 
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NcF 


Colour Factors 
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qa 
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O' 
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- 
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Xi daU 


dg ^ xtK 


1 


1 

2Nc 


- 


- 


10 


1 


K 


H 












11 


1 


d 


H 

Xi UqCI 
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1 
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Table 16: Matrix elements for the Minimal Supersymmetrie Standard Model gaugino 
squark production processes. The number of different colour flows is given by Nqf, ID 
gives the type of diagram, CF gives the colour flow for a given diagram and VP gives the 
virtual particle exchanged in the s-, t- or u-channels depending on the diagram. As before 
the mass eigenstate of the squark is given by a. 



form of the vertex to be ia^P\. The helicity amplitude is given by 



M = ^ , . . (D.l) 
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Figure 34: Two body decay of a fermion to a fermion and a scalar. 
D.1.2 Diagram 2 



Figure 35: Two body decay of an antifermion to an antifermion and a scalar. 



The decay of an antifermion to an antifermion and a scalar is shown in Fig. |35 
In the MSSM this process occurs in the decay of a gaugino to an antifermion and 
a sfermion when we regard the incoming gaugino as an antifermion. The helicity 
amplitude is given by 

M = —= \ ■ [ a~^°s_Ao(/o,Po)i^(Ao,Po,Pi,-"^i,-Ai,/i) (D.2) 

-a^°moF{-\o,lo,Pi,-mi,-\i,li)] . 

As we will usually consider the incoming gaugino to be a particle rather than an 
antiparticle the sign of Aq will be opposite. 

D.1.3 Diagram 3 




Figure 36: Two body decay of a fermion to a fermion and a scalar. 

The decay of a fermion to a fermion and a massless gauge boson, Fig. ^ cannot 
occur in either the Standard or Minimal Supersymmetric Standard Models at tree 
level. However, the processes ~^ X^l and g Xig via loop diagrams can be 
important in the MSSM. 
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The matrix elements for these processes must have the following form [45] 

M = ^n(pi) (Pn - eoe^PL) P2,elu{p,) , (D.3) 

by gauge invariance, where eo is the sign of the mass of the decaying fermion, ei is the 
sign of the mass of the fermion produced in the decay, and (t^*^ = | {'y^'y'^ — t'^t'^)- 
The effective coupling Qes is given in [45] for the decays of the neutralino to a neu- 
tralino and a photon, and in [46] for the decay of a gluino to a neutralino and a 
gluon. 

We can write the helicity amplitude for this decay in terms of the effective 
coupling giving 

M = , ^'"^ (D.4) 

2MoVPo • kPi ■ hP2 ■ k 

[5x2-F{-\i, li,pi, mi, +,p2)s+{p2, l2)F{-,p2,po, mo, -Aq, Iq) 
-5A2+eoeiF(-Ai,/i,pi,mi, -,p2)s-{p2, k2)F{+,p2,po,mo, -Ao,/o)] • 

D.1.4 Diagram 4 




Figure 37: Two body decay of a scalar to a fermion and an antifermion. 

The decay of a scalar to a fermion and an antifermion occurs in the MSSM in 
either the decays of the Higgs bosons or the sfermions. The matrix element is given 

by 

M = - \ == [ a^^Ks^Xjik,Pi)F{\i,pi,p2,-m2,-X2j2) 

2VPi ■ hP2 ■ h 

+a^imiF(-Ai, Zi,p2, -"^2, -A2, ^2)] • 

For the decay of the Higgs bosons to two gauginos we will wish to treat both outgoing 
gauginos as particles in which case the sign of A4 must be changed. 

D.1.5 Diagram 5 

In anomaly mediated SUSY breaking (AMSB) models [47,48] the mass splitting be- 
tween the lightest neutralino and chargino can be very small. The dominant decay 
mode of the lightest chargino is Xi — ^ Xi^r^. The Monte Carlo event generators 
usually contain the three body decay modes X°ud and x~ — ^ X°dii using the 
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Figure 38: Feynman diagram for the decay xf ~^ Xi'^^- 

constituent quark masses. In AMSB models this decay mode, with constituent quark 
masses, is not kinematically accessible which means that the decay of the chargino to 
a pion cannot be described by the three body decay mode followed by the hadroniza- 
tion of the quarks. 

The decay xf ~^ Xi"^^ must therefore be included as a 1 ^ 2 process using 
chiral perturbation theory. The matrix element for this process is 



M = j^uiPi)M^Px^^iPo), (D.5) 



where f^^ is the pion decay constant.^ The helicity amplitude for this process can be 
written as 

M = , (D.6) 

AM^VPo ■ loPi ■ h 

[moa"^° {s^Xi{h,Pi)F{\i,pi,p2, 0, -Aq, Iq) + miF(-Ai, h,p2, 0, -Aq, k)} 
+a^°sxoiPoJo) {s-xAh,Pi)F{Xi,puP2,0,Xo,Po) + miF(-Ai, /i,p2, 0, Ao,Po)}] • 

D.2 Three Body Decay Feynman Diagrams 

There are six Feynman diagrams we will need for the SM and MSSM three body 
decays we consider. Four of these diagrams are for the decays of fermions via either 
virtual gauge boson or scalar exchange, one diagram for the decay of an antifermion 
via gauge boson exchange and one diagram for scalar decay via gauge boson exchange. 
As we are considering all the gauginos to be particles this means we do not need the 
diagrams for antifermion decay via virtual scalar exchange. In all these diagrams Pi 
is the four-momentum of the zth particle, U is the reference vector used to define the 
direction of the ith particle's spin, rrii is the the mass of the ith particle, Aj is the 
helicity of the ith particle and mfj = {pi + pj)'^. 

D.2.1 Diagram 1 

The Feynman diagram for fermion decay via gauge boson exchange is shown in 
Fig.p^. This process occurs for both the top quark, via W exchange, and the MSSM 



^We have defined the pion decay constant in an isospin basis and therefore the pion decay 
constant in [49] should be divided by \/2. 
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3 

Figure 39: Feynman diagram for fermion decay via gauge boson exchange. 



electroweak gauginos, via either W or Z exchange. The form of the first and second 
vertices are ia^'~f^P\ and ib^'j^Px, respectively. It is useful to define a number of 
functions in order to simplify the expression for the amplitude 



= (5aiA2&^'s_Ai(^i,Pi)sai(P2,^2), (D.7a) 

&+- = 5x,-xy'rms^xAp2j2), (D.7b) 

bi\ = -6x,-x,b-^'m2S^xAh,Pi), (D.7c) 

= -5A,A,r^^mim2, (D.7d) 



and 



A 



01 



B 



12 



f^AoAi { a^" ["^o-^-Ai {h,pi)sx^ ipi, lo) - mls^xi ik,Po)sx^ {po, k)] (D.8a) 
+ a~^"momi [s^Xi{k,Po)sx^{poJo) - S-Xi{h,Pi)sxi{piJo)]} 
+ho-Xi { a^°mi [-s_Xi{h,Pi)sxiiPi,Po)s-Xi{PoJo) + mls^XiikJo)] 

+ a"^«mo [s-xAh,Pi)sxiiPi,Po)s-xAPoJo) - mls^xAhJo)] } , 

SxiX2 {b^^ [mls-xi {k,P2)sxr {P2, h) + "^2■S-Al {li,pi)sx, {pi, k)] (D.8b) 

- b'^^mim2 [s-xAh,Pi)sxi{pij2) + S-Ai(/i,P2)sai(P2, ^2)]} 
+5ai-A2 { b^^mi [s-Xiih,Pi)sxiiPi,P2)s-XiiP2j2) + fnls_xiik,k)] 

- b'^^m2 [s^Xiih,Pi)sxi{pi,P2)s-Xi{p2,h) + mls^Xi{k,h)]} 



where Mb and Tb are the mass and width the of the boson, respectively. 
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The amplitude is given by 
M-. 



[D.9) 



l23 

^++ ' 


^ a 




h,Pi, "mi, 


A2,P3j-r { 


^2,P2,P0, 1^0, 






+ a r ( 


-Ai 


li,pi,mi, 


-X2,P2)^ { 


->^2,P3,Po,mo 






{ F( 


-Ai 


h,Pi,mi, 


-A2,P3)F( 


->^2,k,Po,mo, 


~Ao, Iq) 




+ a-^2F( 


-Ai 


h,Pi,mi, 


A2,^2)i^( 


>^2,P3,Po,mo, 


"Aq, Iq)} 




{ a-^^F( 


-Ai 


h,pi,mi, 


A2,/3)i^( 


>^2,P2,Po,mo, 


— Ao, Iq) 






-Ai 


h, Pi, nil, 


-A2,P2)i^( 


->^2,h,Po,mo, 


~Ao, Iq)} 




{ F( 


-Ai 


h,Pi,mi, 


-A2,/3)i^(- 


-A2,/2,P0,"^0, 


~Ao, Iq) 




+ a~^2F( 


-Ai 


li,pi,mi, 


A2,/2)F( 


><2,h,Po,mo, 


— Ao, Iq)} 



^-^01^23 

D.2.2 Diagram 2 




Figure 40: Feynman diagram for fermion decay via Higgs boson exchange. 



The Feynman diagram for fermion decay via Higgs boson exchange is given in 
Fig.^. This process occurs in the MSSM for electroweak gaugino decay into a 
different electroweak gaugino and SM fermions. The form of the first vertex is ia^Px 
and the second vertex ib^P\. It is easiest to express the matrix element for this 
diagram in terms of a function of each of the two vertices 

yXoXi =a^op^^_Xuli,Pi,m,,\o,Po)sx,{po,lo) (D.lOa) 

yx,x, ^ b^'F{-X,,k,p,,m2,Xs,Ps)sx,{P3J3) (D.lOb) 
-b'^''msF{-X2,l2,P2,fn2, -Ag./g). 

The matrix element is given by 

1 1 



M 



AqAi T AA2A3 
V 1 



4VPo ■ kPi ■ hP2 ■ kP3 ■ h - Ml + iT^M^ 
where M$ and r$ are the mass and width of the exchanged scalar, respectively. 



(D.ll) 
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D.2.3 Diagram 3 




Figure 41: Feynman diagram for fermion decay via antisfermion exciiange. 

The Feynman diagram for fermion decay via antisfermion exchange is shown in 
Fig.^. This process occurs in the MSSM for the decay of a gaugino to a different 
gaugino and SM fermions. As before we will take the first vertex to be ia^Px and 
the second vertex to be ib^P\. As this diagram involves scalar exchange it can be 
written in terms of a function for each of the vertices 



yAoA, ^ a^^F{-X2,l2,P2,m2,Xo,Po)sxoiPo,lo) 
+a"^°moF(-A2, h, P2, ?^2, -Ao, ^o), 

-6~^^m3F(-Ai,/i,pi,mi, -X^Js). 



(D.12a) 
(D.12b) 



The matrix element is given by 

1 



M 



1 



T/A0A2T/A1A3 



WPo ■ kPi ■ hP2 ■ I2P3 ■ h mis - + «r$M$ 
where M$ and r$ are the mass and width of the exchanged scalar, respectively. 

D.2.4 Diagram 4 



(D.13) 




Figure 42: Feynman diagram for fermion decay via sfermion exchange. 

The Feynman diagram for fermion decay via sfermion exchange is shown in 
Fig.E^. This process occurs in the decay of a gaugino to SM fermions and a different 
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gaugino. As before we will take the first vertex to be ia^P\ and the second vertex 
to be ib^Px. As this diagram involves scalar exchange it can be written in terms of 
a function for each of the vertices 

^AoA3 ^ a^^F{XoJo,Po,-mo,X3,P3)sxM,h) (D.14a) 

-a"^^m3F(Ao,/o,Po, -nio, -XsJs), 

yx^x, ^ b-^^F{-X2,l2,P2,m2,-XuPi)s-xAPuh) (D.14b) 

-b^^miF{-X2, l2,P2, m2, Ai, h). 

The matrix element is given by 

= ^ ]. Vi''^^V2^^^ , (D.15) 

4a/Po ■ loPi ■ hP2 ■ kP3 ■ h ^12 - Ml + iVq,Mq> i 2 > 

where M<j, and Tq, are the mass and width of the exchanged scalar, respectively. 
D.2.5 Diagram 5 




^3 

Figure 43: Feynman diagram for antifermion decay via gauge boson exchange. 



The Feynman diagram for antifermion decay via gauge boson exchange is shown 
in Fig. This diagram only occurs in the decay of a top antiquark via W exchange 
because we are treating all decaying gauginos as particles. The amplitude for this 
diagram can be expressed using many of the same functions as for diagram one except 
for 

E^^ = 5aoAi { [?^o'5-Ao(^o,Pi)sao(Pi,^i) -"^?-5-Ao(^o,Po)-5ao(Po,^i)] (D.16) 
+a"^imomi [s„Ao(^o,Po)sao(Po, ^1) - S-Ao(^o,Pi)sao(Pi, ^1)]} 
+5ao-Ai { a^'mo [miS__Ao(^o,/i) - s_Ao(^o,Po)sao(Po,Pi)s-Ao(Pi, ^1)] 

-a'^^mi [mls^xoik.h) - S-Ao(^o,Po)sao(Po,Pi)s-Ao(Pi, ^1)] } • 
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The amplitude is given by 
M- 



(D.17) 



(,23 



[ 



F(-Ao, /o,Po, -"^0, A2,P3)^( A2,P2,Pl,-"^l,-Al,/l) 



+ a^' F{-Xo,lo,po, -mo, -\2,P2)Fi->^2,P3,Pi, -"^i, -Ai,/i)] 
+b+-[ ct^' F(-Ao,/o,Po, -"^0, -A2,P3)-^(-A2,/2,Pi, -"^1, -Ai,/i) 



+a ^'^F{-\o,lo,po,-mo, \2,l2)F{ \2,p3,pi, -nii, -XiJi 



+a^' F(-Ao,/o,Po, 



-mo, 
-mo, - 



>^2,h)F{ A2,P2,Pl, -"^1, -Ai,/i) 
-A2,P2)i^(-A2,/3,Pl, -"^1, -Ai,/i 



+623_ [ ir(_^p^/p^pp^_^g^_^2,/3)i^(-A2,/2,Pi 
+a"^"F(-Ao,/o,Po, -"^0, >^2,k)F{ \2,h,Pi: 

+ 1^01^23] 

2 



, -mi, -Ai,/i) 
-mi, -Ai,/i)] 



where and Tb are the mass and width of the gauge boson, respectively. 
D.2.6 Diagram 6 




Figure 44: Feynman diagram for scalar decay via gauge boson exchange. 

The Feynman diagram for scalar decay via gauge boson exchange is shown in 
Fig.^. This process occurs in the decay of both the Higgs bosons of the MSSM, 
to another Higgs boson and SM fermions, and of the sfermions, to another sfermion 
and SM fermions. Again the amplitude for this diagram can be written using many 
of the same functions as for diagram one. The form of the first vertex is ia{po + pi)^ 
and the second vertex ih^^^P\. The amplitude is 



M 



2VP2 ■ hps ■ h "^is + iTbMb 

[ bf^F{X2,P2,P0 + Pl,0, X2,P3) + bf_F{-X2,l2,P0 + Pl,0, -X2,P3 

+b%F{X2,P2,Po + Pi, 0, A2, h) + 6^^_F(-A2, /2,Po + Pi, 0, -A2, /3) 
Hpo+Pi) ■ ij>2+Pz)B^^] ■ 



(D.18) 
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D.3 Four Body Decay Feynman Diagrams 

There is only one four body decay process wliicli we consider. This is the decay of the 
MSSM Higgs bosons to gauge boson pairs, either W or Z, followed by the decay of 
the gauge bosons to fermions. As before pi is the four-momentum of the ith particle, 
li is the reference vector for the ith particle, is the mass of the ith particle. A, is 
the helicity of the ith particle and m^^ = {pi + pj)'^. 




Figure 45: Four body Higgs boson decay to a gauge pair followed by the decay of the 
gauge bosons. 



The Feynman diagram for this process is shown in Fig.^. We will take the 
form of the Higgs boson-gauge boson pair vertex to be iag^'^ , the form of the vertex 
between the first gauge boson and fermion-antifermion to be ih^'^^P\ and the form 
of the vertex between the second gauge boson and the fermion-antifermion to be 
ic^'^f^Px- It is useful to define the functions c^^_j_ and which are the functions Ir"^^ 
and B^"^ with the coupling replaced by c^. 

The helicity amplitude is given by 



M = (D.19) 
a 1 1 

{5aiA3 [ ( SAi(Pi,P3) [4+S_Ai(P4,P2) +c5L + S_Ai(/4,P2)] 

+ ( SAi(P3,P2) [c++S_xAh,Pi) + ci + S_Ai(^l, ^4)] 

+ S-Ai(^l,/3) [c^+-S\^{pi,P2) + C^--Sx^{k,P2)]) 
+ ( SAi(pi,P3) [4+S_Ai(p4, h) + C%S^xAh. h)] 

+ ( Sx,{p3j2) [4+S_Ai(/l,P4) +ci + S-Ai(/l,/4)] 

+ S_Ai(/l,/3) [c+.Sx,{p4j2) + cltsxAh,l2) 
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Decay 


Colour Factor 


Diagram Type 




g ^ qqa 


1 

2 


1 


gq^q 


g qqa 


1 
2 


2 


a^- 
gq^q 


g xh 


1 
4 


3 





Table 17: Two body decays of the gluino. As we are treating the gluino as a particle 
the sign of the gluino's hehcity should be changed in the expression for the second decay 
amplitude. In the decays of the gluino to squarks a gives the mass eigenstate of the squark. 



+ Sx.iPuh) [41s_Ai(P4,P2) + c!is„Ai(^4,P2)]) 
+ ( S^xAh,P3) [c++SxAP4,P2) + C^^+SxAh,P2)] 

+sxAh,P2) [c+^s^xAh,Pi) + clts-x^{k,k)]) 

+ 6i+ ( S_Ai (P3, ^2) [c++Sx^ (Pl, P4) + c!. + SAi ipl, h)] 

+ SxAPl,h) [4-^"Ai(p4,/2) +C^ls_Ai(/4,/2)]) 

+ ( S_Ai(/l,P3) [C++Sx,{p4, I2) + c!.+Sai(/4, ^2)] 

+ SA,(/3,/2) [41s_A,(/l,P4) +c!.ls_A,(/l,/4)])] 

+^512 ^ c:^4^F(A3, P3, Pi + P2, 0, A3, Pa) + 41i^(-A3, /s, Pi + P2, 0, -A3, ^4) 

+c!.^+F(A3,P3,Pl +P2,0,A3,Z4) +c!ii^(-A3,/3,Pl +P2,0,-A3,Z4 
+ [ 5l2^i7^(A,,p,,p3+p4,0,Ai,P2)+&+-i^(-Al,/l,P3 + P4,0,-Ai,P2 

+6i^+F(Ai,pi,p3 + P4, 0, Ai, h) + 6i^_F(-Ai, /i,P3 + P4, 0, -Ai, 

where M^. and F^^ are the mass and width of the zth gauge boson. 
D.4 Matrix Elements 

We can now use the amplitudes for the two and three body decays in order to express 
the matrix elements for decay processes we are interested in. As with the hard 
processes this is simply a matter of listing the Feynman diagrams which contribute 
to a given decay and the various colour factors. This is somewhat easier here because 
the SM and MSSM decays we are considering have only one colour flow so we do not 
have the additional complication of dealing with different colour flows as we did in 
the hard production processes. 

The two body decays are particularly simple to specify because there is only 
one Feynman diagram for each decay mode. The two body decay modes of the 
gluino, neutralinos and charginos are given in Tables p!7| p!8| and IT9|, respectively. The 
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Decay 
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Diagram Type 
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Table 18: Two body decays of the neutralinos. As we are treating the neutralino as a 
particle the sign of the neutralino's helicity should be changed in the expression for the 
second type of decay amplitude. In the decay modes of the neutralino to sfermions a gives 
the mass eigenstate of the sfermion. The colour factor given for sfermion production is 
for the production of sleptons, the corresponding colour factor is Nc for the production of 
squarks. 



two body decays of the Higgs bosons are given in Table pO| and of the sfermions in 
Table. ||. 

The three body decays are more complicated as there is often more than one 
diagram contributing to a given decay mode. The three body SM decay modes of 



the top quarks and antiquarks are given in Table 22. The three body decays of 



the gluinos to quarks and electroweak gauginos are given in Table p3|. The three 
body decays of the neutralinos to fermions and gauginos are given in Table p^. The 



chargino three body decay modes to gauginos are given in Table E5 



The decays of the MSSM Higgs bosons via gauge boson exchange are given 
in Table The decays of the sfermions via gauge boson exchange are given in 
Table |77|. It should be noted that for processes in which the fermions produced in 
these sfermion decays have the same flavour as the decaying sfermion there can be 
additional diagrams involving gaugino exchange which we have neglected. 

We have only considered the four body decays of Higgs bosons, via either real 



or virtual gauge boson pairs. These decay modes are given in Table 28 
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Table 20: Two body decays of the MSSM Higgs bosons. As we are treating the electroweak 
gauginos as particles the sign of the second outgoing particle's helicity should be changed 
for the Higgs boson to two gaugino decay modes. The colour factor given for fermion 
production is for the production of leptons, the corresponding colour factor is Nc for quark 
production. 
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Table 21: Two body decays of the sfermions. Again as we are treating the gauginos as 
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Table 22: Three body top quark decay modes. The colour factors given are for the 
production of leptons in the decays, the colour factor is Nc for the production of quarks. 
The isospin partner of f is f. 
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Table 23: Three body gluino decay modes. As with the hard production cross sections 
the mass eigenstates, a, of the intermediate squarks should be summed over. The isospin 
partner of q is q'. 
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Table 24: Three body neutralino decay modes. We have only included the Higgs boson 
exchange diagrams for neutralino decays. These diagrams are only important for the 
production of bottom and tan due to the Higgs couplings, the production of the top quark 
is is usually kinematically forbidden. In most SUSY models the gluino is heavier than the 
electroweak gauginos and we have therefore not included the modes gqq- As with 

the hard production processes the mass eigenstates of the intermediate sfermions must be 
summed over. The isospin partner of f is f. 
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Table 25: Three body chargino decay modes. The positive chargino can decay to xj^ud 
or ^vl'^ and the negative chargino to y^dn or y^l~v. The gluinos are usually heavier 
than tlic charginos in SUSY models and wc have tlicrcforc not included the decay modes 
xt gqq'- The mass eigenstate, a, of the intermediate sfermions must be summed over. 
As before the isospin partner of a fermion or sfermion is denoted with a prime. 
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Table 26: Three body MSSM Higgs boson decay modes. The colour factors given are for 
the decay of the Higgs boson to a different Higgs boson and leptons. For the decay of the 
Higgs boson to another Higgs boson and quarks the colour factor is Nc- As before f is the 
isospin partner of f. 
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Table 27: Three body decays modes of the sfermions. The fermion n and antifermion 
h produced in the sfermion decay need not be the same as the flavour decaying sfermion. 
The isospin partner of n is denoted by n'. In fact if the produced fermion is of the same 
flavour as the decaying sfermion there are additional Feynman diagrams which we have not 
included. As before the colour factors given are for the production of leptons, the colour 
factor is Nc for the production of quarks. 
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Table 28: Four body Higgs boson decay modes. The colour factor is 1 if all the produced 
fermions are leptons, Nc if a lepton and neutrino and a quark/antiquark pair are produced 
and is all the produced particles are quarks. The isospin partners of n and f are denoted 
by n' and f, respectively. 
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